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I  INTRODUCTION 


The  rapid  proliferation  of  laser  systems  for  civilian  and  military 
applications  urgently  requires  the  development  of  suitable  guidelines  to 
ensure  safety  in  their  use.  In  the  absence  of  appropriate  and  realistic 
guidelines,  excessively  conservative  "safety  factors"  may  be  proposed 
that  could  unduly  restrict  the  use  of  laser  systems. 

The  U.S.  Air  Force  is  actively  conducting  and  sponsoring  research 
studies  to:  (1)  determine  threshold  levels  for  laser  eye  damage 
(Vassiliadis  et  al.f  1969)  and  (2)  determine  the  effects  of  atmospheric 
thermal  turbulence  on  eye-damage  probabilities.  The  latter  is  the  objec¬ 
tive  of  our  present  work  and  the  subject  of  this  report.  Three  previous 
reports  (Johnson  et  al.,  1968,  1970,  and  Dabberdt  and  Johnson,  1971) 
described  the  results  of  the  earlier  research. 

Atmospheric  temperature  inhomogeneities  increase  eye-damage  prob¬ 
abilities  through  the  breakup  of  the  propagating  laser  beam.  To  a 
viewer,  this  appears  as  scintillation  (fluctuations  in  intensity).  The 
"hot  spots,"  or  areas  of  localized  beam  intensification,  typically  have 
intensities  that  are  tens  or  hundreds  of  times  the  average  beam  intensity. 
Thus,  as  the  state  of  the  atmosphere  changes  from  a  thermally  "quiet" 
situation  (no  scintillation)  to  one  of  increasing  thermal  turbulence  and 
scintillation,  the  probability  of  laser  eye  damage  increases  drastically. 

Substantial  progress  has  been  made  toward  clarifying  and  parameter¬ 
izing  the  effects  of  atmospheric  thermal  turbulence  on  laser  propagation 
and  on  the  associated  eye-damage  probabilities.  It  was  determined  exper¬ 
imentally  that  the  amplitude  probability  distribution  of  the  laser  inten- 

*References  are  listed  at  the  end  of  the  report. 
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sity  measured  by  a  5-mm  aperture  detector  is  well  described  by  a  log¬ 
normal  distribution,  in  accordance  with  Tatarski's  (1961)  theory.  The 
results  demonstrated  that  although  the  observed  probability  distributions 
departed  somewhat  from  log  normal  at  the  extreme  upper  end  of  the  dis¬ 
tribution  (possibly  because  of  the  finite  laser  power) ,  the  assumption 
of  a  log-normal  distribution  gave  safe  (conservative)  eye-damage  prob¬ 
ability  estimates.  These  analyses  were  carried  out  to  probabilities 
approaching  10”7,  almost  three  orders  of  magnitude  lower  than  previously 
available . 

With  the  use  of  the  log-normal  model,  a  preliminary  set  of  laser 
safety  guidelines  for  atmospheric  effects  was  developed.  These  guide¬ 
lines  are  based  upon  readily  available,  conventional  weather  observations. 
A  number  of  approximations,  however,  had  to  be  made  in  developing  the 
initial  guidelines,  because  of  the  basic  lack  of  information.  Additional 
experiments  were  conducted  to  clarify  certain  unresolved  propagation 
effects  and  to  extend  the  guidelines  to  a  broader  range  of  laser  system 
use  in  terms  of  wavelength,  path  geometry,  nature  of  the  platform,  and 
turbulence  conditions.  Our  results  indicated  the  presence  of  a  super¬ 
saturation  region  of  optical  scintillations  along  a  horizontal  path  in 
which  the  scintillation  magnitude  decreases  with  increasing  range  and 
thermal  turbulence  level.  The  eye-safety  guidelines  were  revised  to 
incorporate  this  effect  as  well  as  the  results  of  our  initial  scin¬ 
tillation  measurements  along  a  slant  path.  Additional  experimental  work 
was  carried  out  to  investigate  slant-path  propagation  and  the  wavelength 
dependence  of  scintillation.  The  results  of  this  research  were  also 
incorporated  into  the  eye-safety  guidelines. 

During  the  past  year,  further  experimental  work  was  carried  out  to 
investigate  in  greater  detail  the  nature  of  the  scintillation  for  a  slant 
propagation  path.  Measurements  were  made  simultaneously  along  identical 
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ground-to-air  and  air-to-ground  paths.  Concurrent  vertical  profiles  of 

thermal  turbulence  (i.e.,  the  refractive  index  structure  function,  C  ) 

n 

were  also  obtained  and  have  been  used  to  parameterize  atmospheric 
effects.  The  results  of  the  study  indicate  the  onset  of  scintillation 
saturation  at  long  ranges,  low  elevation  angles,  and  high  thermal  tur¬ 
bulence.  A  significant  finding  of  our  experimental  work  has  been  the 
heretofore  unobserved  dependence  of  the  scintillation  pattern  on  the 
scale  of  atmospheric  thermal  turbulence,  in  addition  to  the  widely  recog¬ 
nized  dependence  on  range,  path  geometry,  and  the  vertical  profile  of  the 
refractive  index  structure  function.  More  specifically,  when  the  slant- 
path  scintillation  (at)  predicted  by  the  classical  theory  of  Tatarski 
(1961  and  1971)  is  taken  us  the  scaling  factor  for  range,  geometry,  and 

the  vertical  profile  of  C  ,  the  observed  scintillation  shows  evidence  of 

n 

saturation  at  o^-values  that  are  about  a  factor  uf  two  larger  during 
unstable  (lapse)  atmospheric  conditions  than  during  stable  (inversion) 
conditions.  The  pattern  is  similar  for  both  paths,  although  the  observed 
scintillation  at  saturation  is  slightly  larger  (about  1.3  dB)  for  the 
ground-to-air  path.  These  results  and  the  incorporation  of  this  infor¬ 
mation  into  the  eye-safety  guidelines  are  described  in  the  main  body  of 
this  report. 

Additionally,  we  have  developed  a  simplified  eye-safety  evaluation 
procedure  for  the  estimation  of  the  worst-case  probability  of  eye  damage. 
The  worst-case  estimate  is  made  using  only  the  laser  and  operational 
characteristics,  while  the  atmospheric  effects  are  treated  implicitly  in 
the  procedure.  This  worst-case  guide  is  intended  for  use  when  operational 
restrictions  prohibit  the  use  of  the  more  detailed,  general  guidelines. 
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II  PROGRAM  SUMMARY 


A.  Background 

Classical  optical  propagation  theory  (Tatarski,  1961  and  1971) 

predicts  that  the  scintillation  magnitude  for  a  near-ground,  horizontal 

propagation  path  is  directly  proportional  to  the  refractive  index 

structure  function  (C  ) ,  approximately  directly  proportional  to  the 

n 

range  (R) ,  and  approximately  inversely  proportional  to  the  square  root 
of  the  optical  wavelength  (X);  the  exact  form  is 

(a  )  =  K  R11/12  (2tt/X)7/12  C  .  (1) 

t  H  1  n 

Here  a  is  the  theoretical  standard  deviation  of  the  log  signal  intensity, 

X 

K,  is  a  dimensionless  constant ,  and  the  subscript  H  denotes  a  horizontal 
path.  This  theory,  which  is  based  on  the  so-called  Rytov  approximation, 
fits  observations  well  for  low  thermal-turbulence  conditions  and  rela¬ 
tively  short  ranges. 

Tatarski 's  theory  also  predicts  the  scintillation  magnitude  for  the 
case  of  optical  propagation  along  a  slant  path  (subscript  S) , 

r  ^  11/2 

(ot)s  =  \  (2n/\>7/12  f  c2  (r)  rV6  dr  J  ,  (2) 

where  the  coordinate  origin  (r  =  0)  is  at  the  receiver.  The  term  shown 

in  brackets  is  referred  to  as  the  range  integral  and  is  a  function  of 

range,  path  geometry,  the  slant  profile  of  thermal -turbulence  intensity, 

and  receiver-transmitter  orientation.  For  the  normal  case  of  a  decrease 

in  C  with  increasing  height,  the  theory  predicts  that  the  scintillation 
n 
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will  be  larger  for  a  ground-to-air  propagation  path  than  vice  versa.  The 
results  of  our  earlier  slant-path  studies  (Johnson  et  al.,  1970,  and 
Dabberdt  and  Johnson,  1971)  indicated  that  the  applicability  of  Equation 
(2)  may  be  restricted  to  relatively  small  values  of  the  range  integral. 

Application  of  the  theory  to  longer  ranges  and  higher  turbulence 
levels  has  not  been  Intuitively  appealing  since  for  these  conditions,  the 
predicted  scintillation  magnitude  (and  hence  the  eye-damage  probability) 
increases  indefinitely  without  bound.  However,  no  better  theory  has  been 
available,  and  indeed,  recent  advances  have  come  principally  through 
experimentat ion . 

The  first  breakthrough  came  from  experiments  in  Russia  by  Gracheva 
and  Gurvich  (1965)  and  by  Gracheva  (1967),  which  showed  a  saturation,  or 
leveling  off,  of  the  scintillation  magnitude  with  increasing  range  beyond 
about  1  km.  (The  longest  range  at  which  measurements  are  reported  in  the 
Russian  work  is  1.75  km.)  Subsequently,  Tatarski  (1966)  and  DeWolf  (1968) 
developed  variations  of  the  propagation  theory  that  conform  to  some  degree 
to  the  Russian  observations.  A  second  breakthrough  resulted  from  our 
earlier  work  (Johnson  et  al.,  1970)  at  a  field  site  in  Woodland,  California, 
that  showed  the  presence  of  a  supersaturation  region  where  the  scintil¬ 
lation  magnitude  for  a  helium-neon  laser  and  a  near-ground,  horizontal  path 
actually  decreases  with  increasing  range  and  thermal-turbulence  intensity 
beyond  the  saturation  level.  An  analytical  expression  was  given  that  relates 

the  measured  scintillation  magnitude  (o  )  to  the  theoretical  (a  )  : 

m  H  t  H 

<at’H 

(a  >  =  1  +  0.16  (3) 

m  H  t  H 


Figure  1  gives  a  comparison  of  the  measured  values  of  the  scintillation 
intensity  with  the  values  predicted  by  Tatarski 's  theory. 
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FIGURE  1  COMPARISON  OF  MEASURED  VALUES  OF  LOG  INTENSITY 
STANDARD  DEVIATION  WITH  VALUES  PREDICTED  BY 
TATARSKI'S  THEORY,  SHOWING  THE  SUPERSATURATION 
EFFECT,  FOR  THE  WOODLAND  EXPERIMENT  (Johnson  et  «l„ 
1970) 


Another  experimental  program  was  conducted  in  September  1970,  at  an 
excellent  field  site  in  Mercury,  Nevada,  that  was  designed  to  investigate 

the  dependence  of  optical  scintillations  on  range,  thermal-turbulence 

intensity,  path  geometry,  and  wavelength.  Because  of  the  location  and 

time  of  the  experiment,  a  significantly  wider  variation  (~  50%)  of  C^ 

values  were  measured  over  longer  ranges  (to  7.6  km)  than  in  our  earlier 

experiments.  We  found  that  the  measured  scintillation  magnitude  for  a 

near-ground  horizontal  path:  (1)  saturated  with  Increasing  range  and 

thermal-turbulence  level  (as  represented  by  the  scaling  factor  o.„) , 

vH 

(2)  reached  this  peak  (9  *)  at  about  9.u  ■  2.5  and  decreased  thereafter 

Bln  tH 

("supersaturation"),  and  (3)  showed  evidence  of  leveling  off  beyond 
a  m  IS  at  a  value  of  approximately  0.5.  An  analytical  expression  that 

vH 

fits  the  0.6328-m.  data  well  is 
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(4) 


(Vh 


O  (0.6328  p)  =  1  +  0.30  (a  )  5/3 
mH  t  H 


where  a  is  the  measured  log-intensity  standard  deviation  along  a 
mH 

horizontal  path.  Furthermore,  our  multiwavelength  measurements  (0.4880  p, 

* 

0.6328  p,  and  1.064  p)  showed  that:  (1)  the  ratios  of  a  for  the  various 

mH 

wavelengths  are  well  approximated  by  Tatarski’s  "minus  seven-twelfths" 

*■ 

dependence  on  wavelength;  (2)  the  value  of  a  at  the  measured  peak  (o  ) 

tH  mH 

varies  inversely  with  wavelength;  and,  most  significantly,  (3)  optical 
scintillations  show  little  apparent  dependence  on  wavelength  in  the  far- 
supersaturation  region.  Analytical  expressions  for  the  0.4880-p  and 
1.064-p  data  are  given  by 


and 


a  (0.4880  p)  =  1  +  0.25  (.a  )  1,77 
mH  t  H 


amH  (1,064  u)  =  1  +  0.48  (  7  )  1,60 

t  H 


(5) 


(6) 


Application  of  Tatarski's  slant-path  theory  is  also  restricted  owing 
to  the  monotonic  increase  of  the  predicted  scintillation  magnitude  with 
increasing  turbulence  levels  and  range.  The  results  of  the  earlier  pro¬ 
gram  (Johnson  et  al.,  1970)  indicated  the  possibility  of  scintillation 
saturation  along  a  slant  path,  although  the  limited  range  (^  4000  m)  pre¬ 
cluded  possible  supersaturation.  During  the  1970  field  study,  static  slant- 
path  scintillation  measurements  were  made  using  a  465-m  tower  as  a  receiver 
platform  for  a  ground-based  helium-neon  laser  (0.6328  p) ;  vertical  profiles 
of  thermal  turbulence  intensity  were  also  collected.  These  measurements 
were  made  during  daytime— lapse  conditions— and  the  measured  scintillation 
magnitude  (a  )  was  observed  to:  (1)  saturate  with  increasing  range  and 

thermal-turbulence  level  (as  represented  by  a  )  and  (2)  show  some  indi- 

to 

cation  of  supersaturation  beyond  o  a*  8.  The  observed  scintillation 

ts 
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maxima  for  the  slant-  and  horizontal-propagation  paths  were  virtually 
the  same  at  about  1.25.  We  found  that,  for  our  experimental  conditions, 
the  eye-damage  hazard  was  larger  along  the  slant  path  for  ranges  greater 
than  4  km.  This  resulted  from  saturation  along  the  slant  path  together 
with  supersaturation  on  the  horizontal  path  at  these  ranges. 

B.  Summary  of  Results 

An  extensive  field  program  was  conducted  during  July  197.1  at  the 
Mercury  site  to  study  the  nature  of  the  scintillation  of  a  helium-neon 
laser  beam  for  ground-to-air  and  air-to-ground  propagation  paths  during 
both  daytime  (lapse)  and  nighttime  (inversion)  periods.  The  scintillation 
measurements  were  made  simultaneously  along  the  two  coincident  paths, 
together  with  concurrent  measurements  of  the  differential-temperature 

structure  (used  to  determine  the  refractive  index  structure  function,  C  ) 

n 

at  four  heights  from  the  surface  to  460  m  (the  height  of  the  elevated  laser 
and  receiver).  Supplemental  mean  wind  and  temperature  profiles  were  also 
obtained. 

One  significant  aspect  of  the  program  was  the  data  it  provided  on  the 

time-height  variations  of  the  refractive  index  structure  function.  As 

expectec,  C  is  a  minimum  during  the  near-sunrise  and  -sunset  adiabatic 
n 

periods,  while  nighttime  (inversion)  values  are  typically  about  one-half 

those  observed  during  the  midday  (lapse)  periods.  The  general  feature  of 

the  vertical  C  profile  is  the  quaBi-exponential  decay  with  height.  How- 
n 

ever,  during  the  sunrise  and  sunset  stauillty-transition  periods,  a  secondary 
Cn  maximum  occurs  well  above  the  surface.  This  feature  is  particularly  sig¬ 
nificant  when  evaluating  the  slant-path  scintillation  in  terms  of  the  theo¬ 
retical  value  given  by  Equation  (2).  A  simple  model  that  describes  the 
features  of  the  time-height  dependence  of  is: 

C  (z,t)  =T'{c  +  iC  cos  (w  t  -  z/Zt  +  £)}  ,  (7) 

n  n,i  in  i  w  * 
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,y,  i 
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The  overbar  denotes  the  time  average  of  the  i  .  harmonic,  AC  is  the 

th  n 

harmonic  amplitude,  uj  is  the  frequency,  Z  is  the  damping  depth,  and  £  is 
the  phase  lag. 

A  particularly  significant  finding  of  the  program  is  the  apparent 
stability  dependence  of  the  scintillation  pattern.  As  with  the  earlier 
horizontal  path  studies,  a  [in  this  case,  Equation  (2)]  is  used  as  a 

T* 

range-turbulence-wavelength  scaling  parameter.  While  the  observed  scin¬ 
tillation  patterns  for  both  the  ground-to-air  and  air-to-ground  paths  are 
similar  for  the  two  stability  categories  (see  Figures  2  and  3),  corres¬ 
ponding  absolute  values  of  a  are  twice  as  large  for  the  lapse  cases  than 
for  the  inversion  cases.  This  suggests  that  there  is  a  scaling  effect 
not  only  of  range,  geometry,  wavelength,  and  turbulence  intensity,  but 
also  of  the  stability  types,  i.e.,  lapse  or  inversion.  We  therefore  pro¬ 
pose,  in  way  of  conjecture,  that  the  scintillation  magnitude  is  dependent 
on  the  frequency  (or  scale)  as  well  as  the  magnitude  of  the  thermal 
fluctuations. 

During  lapse  conditions,  there  is  a  decrease  in  the  atmospheric 
potential  temperature  with  height  and  a  tendency  for  convective  activity 
resulting  from  the  unstable  density  stratification;  for  inversion  conditions 
the  reverse  is  true  and  the  atmosphere  is  regarded  as  stable.  Typically, 
the  spectrum  of  near-surface  thermal  fluctuations  during  inversion  con¬ 
ditions  (typical  of  clear  nights)  is  characterized  by  relatively  high  fre¬ 
quencies  (and  low  wavelengths)  due  to  the  height  variation  of  temperature 
and  the  action  of  mechanical  turbulence.  During  lapse  conditions  (char¬ 
acteristic  of  sunny  days),  the  spectrum  of  thermal  fluctuations  is  often 
augmented  by  the  effects  of  buoyancy  (free  convection).  The  impact  of  this 
convective  activity  is  reflected  in  the  occurrence  of  a  second  peak  in  the 
spectrum  at  lower  frequencies  (higher  wavelengths).  In  summary,  nighttime 
conditions  are  characterized  by  relatively  small  though  frequent  turbulent 
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FIGURE  2  COMPARISON  OF  MEASURED  VALUES  OF  LOG 
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PREDICTED  BY  TATARSKI'S  THEORY,  SHOWING 
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elements,  whereas  daytime  conditions  with  free  convection  may  be  dominated 
by  the  larger,  infrequent  convective  cells.  As  a  consequence,  it  may  be 
that  the  "effective  range"  (or  number  of  refructive  occurrences)  is  sig¬ 
nificantly  less  for  these  daytime  conditions  although  the  value  of  C  may 

n 

not  differ  from  the  nighttime  case.  The  idea  of  an  effective  range  is 

introduced  to  specify  the  relationship  between  the  scale  of  turbulence 

and  the  scintillation:  For  equal  values  of  C.  and  a  given  propagation 

n 

path,  daytime  scintillation  under  free  convection  is  observed  to  occur  as 
if  the  range  were  effectively  reduced. 

The  occurrence  of  scintillation  saturation  is  indicated  at  large  a 
for  all  path  and  stability  categories,  but  there  is  no  evidence  of  a 
supersaturation  effect  within  the  range  of  a  .  values  observed.  Again 

there  appears  to  be  an  influence  of  the  stability  category— -in  this  case 

* 

on  the  magnitude  of  the  scintillation  at  saturation,  a  •  For  the  irver- 
*  m 
sion  cases,  o  is  about  1.7,  while  for  the  lapse  cases  it  is  about  1.4. 
m 

This  difference,  though  only  about  20%  or  1.3  dB,  is  quite  important  in 
the  evaluation  of  the  eye-damage  probability  because  of  non-linear 
effects. 

I 

These  findings  have  been  incorporated  into  a  revised  procedure  for 
the  evaluation  of  the  probability  of  laser  eye-damage.  Additionally,  a 
short  and  simplified  procedure  has  been  developed  to  estimate  the  worst- 
case  probability  of  eye  damage. 


\ 
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Ill  EXPERIMENTAL  PROGRAM 


A .  General 

The  general  objectives  of  the  program  were  to  extend  and  supplement 
our  previous  laser  eye-safety  research  (Johnson,  Evans,  and  Uthe,  1970, 
and  Dabberdt  and  Johnson,  1971)  and  to  develop  a  simplified  laser  eye- 
safety  guide  covering  a  broad  range  of  use  of  laser  systems  in  terms  of 
path  geometry,  atmospheric  conditions,  and  laser  characteristics.  Specif¬ 
ically,  the  objectives,  of  the  experimental  program  were  to: 

(1)  Verify  the  existence  of  scintillation  saturation  and 
supersaturation  for  ground-to-air  propagation,  and 
simultaneously  determine  whether  saturation  and 
supersaturation  occur  for  the  air-to-ground  path. 

(2)  Examine  the  path  dependence  of  the  magnitude  of  the 
scintillation  maximum. 

(3)  Develop  relationships  between  slant-path  scintillation 
and  atmosphoric  structure. 

Toward  these  ends,  an  experimental  program  was  conducted  during  the 
period  8  to  23  July  1971  at  the  Nevada  Test  Site  (NTS).  Simultaneous 
measurements  were  made  sequentially  of  the  intensity  fluctuations  in  the 
beams  of  two  15-mW  helium-neon  lasers,  along  ground-to-air  and  air-to- 
ground  propagation  paths  at  six  ranges.  A  chronological  summary  of  our 
activities  is  given  in  Table  1.  Two  receivers  with  5-mm  diam.  apertures 
were  used,  one  near  the  top  of  the  405-m  BREN  Tower  and  the  >ther  in  a 
mobile  instrumentation  laboratory.  The  horizontal  range  between  the 
tower  and  the  six  sites  varied  between  0.625  and  10.0  km;  both  daytime 
and  nighttime' measurements  were  made.  Coincident  thermal-turbulence 
measurements  were  made  at  four  heights  ranging  from  1.7  to  460  m; 
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Table  1 


Date 

6  July 

7  July 

8  July 

9  July 

10  July 

11  July 

12  July 

13  July 

14  July 

15  July 

16  July 
17-18  July 

19  July 


CHRONOLOGICAL  SUMMARY,  1971  BREN  TOWER 
EXPERIMENT,  NEVADA  TEST  SITE  (NTS) 


Activity 

Equipment  departn  SRI 

Equipment  and  field  party  am  ve  Lae  Vegas 

Equipment  and  field  par+v  arrive  NTS,  clear 
security;  unpack  equipment ,  locate  field 
sites,  calibrate  and  begin  installation  of 
equipment  with  aid  of  EG&G  personnel;  Air 
Resources  Laboratory  personnel  activate  wind 
and  temperature  Bystem  on  BREN  Tower 

Continue  equipment  check-out,  calibration, 
and  installation 

Complete  equipment  installation  and  11  daytime 
data  runs 

Rest  day 

Complete  26  daytime  data  runs 

Begin  nighttime  data  collection  sequence;  13 
runs  completed 

Complete  nighttime  data  collection  sequence; 

19  runs  completed 

Complete  28  daytime  data  runs 

Complete  15  daytime  data  runs 

Rest  days 

Begin  nighttime  data  collection  sequence; 

18  runs  completed 
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Table  1  (Cont.) 


Da  t  e 


20  July 

Complete  nighttime  data  oolleotion 
17  runs  completed 

sequence ; 

21  July 

Complete  9  daytime  data  runs 

22  July 

Complete  26  daytime  data  runs; 

shut 

down 

23  July 

Begin  equipment  take-down  and 

packaging 

24  July 

Complete  equipment  packaging, 
security  and  depart 

clear 

NTS 

26  July 

Equipment  arrives  SRI 

supplemental  wind  and  temperature  profiles  ware  also  made  using  available 
tower  instrumentation. 

B.  Site  Description 

The  BREN  Tower  facility  at  the  NTS  is  looated  about  80  miles  north¬ 
west  of  Las  Vegas  at  Lat  38°  46'  50.4"  N,  Long  116°  14'  33.4"  W,  at  an 
elevation  of  1110  m  msl.  The  tower  is  an  equi triangular  structure  approx¬ 
imately  4  m  on  a  side,  with  an  open  steel  framework.  The  terrain  slope 
within  a  2.5-km  radius  of  the  tower  is  approximately  1:40,  with  the  fall 
line  ENE  to  SSW,  except  for  a  ridge  5  km  distant  in  the  SE  quadrant  (see 
Figure  4) .  All  runs  were  made  over  paths  to  the  west  of  the  tower  to 
ensure  maximum  homogeneity  of  the  surface,  thereby  permitting  us  to  be 
reasonably  confident  that  the  meteorological  conditions  measured  at  a 
fixed  point  over  suitable  averaging  time  would  be  representative  of  the 
entire  propagation  path.  Figure  5  is  a  topographic  cross  section  between 
the  BREN  Tower  and  Site  6  (range  10  km);  the  slope  (1:62)  is  quite  uniform 
over  the  entire  path.  The  surface  is  basically  a  sandy,  pebbly  gravel  of 
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FIGURE  4  MAP  OF  THE  FIELD  SITE  AT  THE  NEVADA  TEST  SITE 


elevation 


volcanic  rock  fragments,  primarily  tuff  with  some  basalt  and  andesite. 

This  uppermost  stratum  (caliche)  extends  to  a  depth  on  the  order  of  30  m. 
The  instrumentation  runway  on  which  the  shorter-range  measurements  (0.63 
km)  were  made  is  devoid  of  vegetation;  the  surrounding  terrain  has  a 
typically  sparse  growth  of  sagebrush. 

C.  Instrumentat ion 

1.  Laser  Measurements 

The  tower  and  truck  laser-receiver  systems  were  virtually 
identical  (Figures  6  and  7).  Two  Spectra-Physics  Model  124  helium/ 
neon  lasers  (wavelength  0.6328  p,)  with  an  output  power  of  approximately 
15  mW  were  used  for  this  experiment;  the  surface  laser  was  positioned 
1.7  m  above  the  ground,  while  the  elevated  laser  was  at  the  460-m  level 
on  the  tower.  Transmitter-receiver  alignment  was  controlled  with  micro¬ 
positioners  on  both  the  laser  mount  and  the  receiver  mount;  a  corner  cube 
retroref lector  mounted  on  the  receiver  aided  in  positioning  the  laser. 
Beam-expanding  and  collimating  optics  (S-P  Model  336)  were  used  to  obtain 
a  transmitted  beam  diameter  of  approximately  20  mm  (measured  between  the 
e"2  intensity  levels).  A  spatial  filter  (S-P  Model  332)  with  an  aperture 
of  15  p,  (S-P  Model  A6)  was  used  to  ensure  a  smooth,  mean  (Gaussian) 
intensity  profile  across  the  beam. 

Two  virtually  identical  receivers  were  used  to  measure  the  laser 
intensity  fluctuations.  One  receiver  was  positioned  near  the  truck  to 
monitor  the  signal  from  the  tower-mounted  laser,  while  the  second  receiver 
was  at  the  460-m  level  on  the  tower  and  monitored  the  signal  from  the  truck 
laser.  Each  receiver  employed  a  5-mm  diam.  aperture  to  simulate  the  average 
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size  of  the  pupil  of  the  human  eye,  a  10-jJ  interference  filter  centered 
at  0.6328  \i  to  reduce  background  (ambient)  radiation,  and  a  3-mrad  full- 
angle  field  of  view.  An  RCA  Model  7265  photomultiplier  tube  (PMT)  with 
an  S-20  photocathode  surface  was  used  in  each  receiver. 

The  output  signal  from  each  PMT  was  fed  into  a  logarithmic 

t 

amplifier  (PM  Model  1002)  and  the  output  voltage  was  then  recorded  in 
FM  mode  on  an  analog  instrumentation  tape  recorder  (H-P  Model  3960A). 

The  mean  PMT  current  was  monitored  by  a  meter  using  an  R-C  circuit  with 
a  time  constant  of  approximately  7  seconds.  Data  from  the  truck-based 
PMT  were  recorded  in  the  truck;  power  for  the  mobile  instrumentation  was 
supplied  by  a  bank  of  storage  batteries.  These  were  capable  of  8  to  12 
hours  continous  operation  and  were  recharged  during  off-periods.  Data 
from  the  tower-based  PMT  were  recorded  in  a  bunker  (see  Figure  8)  at  the 
base  of  the  tower. 

2.  Meteorological  Measurements 

The  primary  meteorological  instrumentation  consisted  of  four 
fast-response  differential  thermometers  (Figure  6)  for  the  measurement  of 
the  temperature  structure  function.  A  vane  was  used  to  maintain  the 
orientation  of  two  horizontally  separated  temperature  sensors  perpendicular 
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to  the  wind  direction;  two  0.0005-in.  diam.  chromel-constantan  thermo¬ 
couples  were  used,  physically  separated  by  60  cm  but  electrically  connected 
in  series.  The  voltage  generated  by  the  temperature  difference  between  the 
thermocouples  (64  u,V/°C,  as  determined  from  repeated  calibration  tests)  was 
then  amplified  by  a  high-gain  differential  amplifier  (designed  and  con¬ 
structed  at  the  Institute)  and  recorded  in  FM  mode  on  the  analog  magnetic 
tape  recorders  (H-P  Model  3960A).  The  time  constant  of  the  differential 
thermometer  was  determined  experimentally  to  be  0.015  s  when  ventilated  at 
2  m/s.  This  response  was  sufficiently  fast  so  that  the  error  in  the  com¬ 
puted  fluctuation  variance  caused  by  neglecting  the  higher  frequencies 
should  be  small. 

The  four  differential  thermometers  were  used  to  provide  a  ver¬ 
tical  profile  of  the  temperature  structure  function  for  the  slant-path 
scintillation  measurements.  The  lower-most  sensor  was  positioned  near  the 
truck  atop  a  tripod  at  the  approximate  height  of  the  surface  laser  (1,7  m) ; 
the  other  three  sensors  were  mounted  on  booms  that  extended  3  m  outward 
from  the  tower  at  heights  of  91.5,  213,  and  460  m.  The  near-surface 
differential  thermometer  was  used  at  the  location  of  the  near-surface  laser 
to  provide  thermal-turbulence  measurements  at  the  point  where  turbulence 
effects  on  the  intensity  fluctuations  are  most  pronounced — near  the  source. 
For  the  elevated  laser,  the  460-m  differential  thermomoter  provided  the 
resolution  required.  Thermal-turbulence  measurements  at  the  two  inter¬ 
mediate  levels  were  used  to  help  defino  the  vertical  structure  of  the 
temperature  structure  function. 

Supplementary  mean-temperature  profiles  were  obtained  as  hard¬ 
copy  output  from  an  electromechanical  printer  at  10-mln.  intervals  from 
existing  National  Weather  Service  (NWS)  sensors  mounted  at  ten  levels  on 
the  tower;  3,0,  16,4,  25,9,  47.3,  90,0,  133,  175,  213,  303,  and  460  m. 

Mean  wind  speed  and  direction  were  taken  visually  from  strip  chart  records 
over  the  course  of  each  run;  NWS  wind  sensors  at  eight  levels  were  used: 
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3.0,  16.4,  25.9,  90.0,  175,  218,  303,  and  460  m.  Additional  measurements 
of  wind  speed  and  direction  were  obtained  using  a  propeller-vane  maintained 
at  a  height  of  1.7  m  at  a  location  approximately  100  m  SW  of  the  tower 
(see  Figure  9).  These  data  were  recorded  on  an  incremental,  dual-channel 
chart  recorder. 

D.  Experimental  Procedure 

The  basic  configuration  of  the  path  geometry  and  range  dependence 
experiment  is  illustrated  in  Figure  10.  Lasers  and  receivers  were  located 
at  the  top  of  the  BREN  Tower  and  sequentially  at  various  horizontal  ranges: 
0.63,  1.4,  2.5,  5.0,  7.5,  and  10.0  km.  The  signal  inputs  from  the  tower 
receiver  and  three  differential  thermometers  were  brought  into  a  recording 
bunker  located  near  the  base  of  the  tower.  A  second  recording  system  in 
the  truck  was  used  for  the  ground-based  laser  scintillation  and  temperature 
fluctuation  measurements. 

The  experimental  procedure  entailed  the  simultaneous  measurement  of 
laser  scintillations  over  reciprocal  slant -propagation  paths  sequentially 
at  the  six  different  ranges;  temperature  fluctuation  data,  together  with 
the  wind  and  temperature  profiles,  were  recorded  simultaneously  with  the 
laser  data.  In  practice,  the  technique  entailed  the  transport  of  the 
mobile  instrumentation  system  to  the  desired  range.  A  6-min.  run  was  made 
at  a  site,  including  calibration  levels  and  a  check  of  the  background 
(ambient  light)  noise.  The  system  was  then  transported  by  truck  to  the 
nexv  site,  and  the  procedure  was  repeated;  transit  and  setup  time  varied 
typically  from  15  to  35  min.,  depending  on  the  range  change  and  the  laser- 
receiver  alignment  time. 

Before  each  run,  the  PMT  supply  voltage  was  adjusted  to  give  an 
approximate  mean  output  current  of  10  uA.  The  log  amplifiers  were  adjusted 
to  give  zero  output  voltage  from  an  input  of  10  pA,  so  that  the  means  of 

the  recorded  (logged)  signals  would  be  near  zero. 

22 


FIGURE  8  RECORDING  INSTRUMENTATION  FOR  TOWER 

RECEIVER  AND  DIFFERENTIAL  THERMOMETERS: 

(1)  FM  INSTRUMENTATION  TAPE  RECORDER, 

(2)  PROBABILITY  DENSITY  ANALYZER,  AND 

(3)  X-Y  PLOTTER 


FIGURE  B  VIEW  TO  THE  WEST  OF  BREN  TOWER  ILLUSTRATING 
THE  TOPOGRAPHY  AND  SHOWING  THE  LOW-LEVEL 
WIND  SENSOR 


FIGURE  10  SCHEMATIC  OF  EQUIPMENT  INSTALLATION 


The  high~gain  differential  amplifier  for  the  near-surface  differ- 

o 

ential  thermometer  was  set  so  that  a  50-p,V  input  signal  (~  0.78  C  differ¬ 
ential)  produced  a  1-V  output.  The  other  three  amplifiers  produced  a 
1-V  dc  output  for  a  25-p.V  input  signal  (~  0.39°C  differential)  owing' to 
the  smaller  temperature  fluctuations  at  the  higher  levels. 

A  total  of  182  data  runs  were  collected  at  six  ranges  for  all  hours 
of  the  day.  The  runs  are  tabulated  in  Table  2. 

E.  Preliminary  On-Site  Data  Analysis 

As  mentioned  earlier,  the  laser  scintillation  and  temperature 
fluctuation  data  were  recorded  in  analog  fashion  on  magnetic  tape.  In 
addition,  the  signals  from  the  tower  instrumentation  were  monitored  on 
line  with  a  probability  density  analyzer  (PDA,  see  Figure  8)  to  examine 
the  quality  of  the  data  in  real  time.  The  PDA  automatically  "builds  up" 
a  probability  density  histogram  (PDH)  as  data  are  being  collected.  When 
a  preselected  number  of  samples  has  been  collected,  the  PDH  is  auto¬ 
matically  displayed  and  can  then  be  recorded  on  an  on-line  x-y  plotter. 
(This  procedure  has  since  been  Improved,  and  the  PDH  is  now  automatically 
transferred  to  a  desk  computer  for  statistical  analysis— see  Chapter  IV). 
In  this  manner  we  were  able  to  monitor  laser-receiver  alignment,  back¬ 
ground  noise  level,  thermal  turbulence,  and  so  forth.  Then,  for  example, 
if  the  mean  laser  signal  level  varied  during  a  run,  we  were  able  to 
remedy  ihe  problem  and  repeat  the  run. 
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Table  2 


DATA  SUMMARY— RECIPROCAL  SLANT  PATH  EXPERIMENT, 
NEVADA  TEST  SITE,  1971 


Site 

ranges: 

1  =  0. 

63  km, 

2  =  1.4  km, 

3  =•  2. 

5  km, 

4=5. 

0  km, 

5  =  7,5  km, 

6  =  10. 

0  km 

RUN 

TIME 

SITE 

RUN 

TIME 

SITE 

NO. 

DATE 

(PDT) 

NO. 

NO. 

DATE 

(PDT) 

NO. 

A 

7/10 

1011 

1 

16 

7/12 

1737 

4 

B 

1025 

1 

17 

1756 

5 

C 

1034 

1 

18 

1812 

6 

D 

1111 

2 

19 

1833 

4 

E 

1132 

3 

20 

1854 

4 

F 

1205 

4 

21 

1915 

4 

G 

1213 

4 

22 

1933 

4 

H 

1238 

6 

23 

7/13 

1910 

1 

I 

1255 

6 

24 

1944 

4 

J 

1325 

5 

25 

1955 

4 

K 

1333 

5 

26 

2012 

4 

L 

7/12 

1012 

4 

27 

2030 

4 

M 

1028 

4 

28 

2057 

4 

N 

1059 

6 

29 

2132 

6 

0 

1142 

6 

30 

2152 

5 

1 

1155 

6 

31 

2211 

4 

2 

1238 

1 

32 

2233 

3 

3 

1247 

1 

33 

2243 

3 

4 

1312 

1 

34 

2301 

2 

3 

1340 

2 

35 

2331 

1 

6 

1406 

3 

36 

7/14 

0026 

6 

7 

1438 

4 

37 

0033 

6 

8 

1446 

4 

38 

0050 

6 

9 

1509 

5 

39 

0109 

5 

10 

1518 

5 

40 

0131 

4 

11 

1538 

6 

41 

0149 

3 

12 

1547 

6 

42 

0205 

2 

13 

1626 

1 

43 

0227 

1 

14 

1653 

2 

44 

0302 

6 

15 

1712 

3 

45 

0327 

5 
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Table  2  (Cont.) 


DATA  SUMMARY- -RECIPROCAL  SLANT  PATH  EXPERIMENT, 
NEVADA  TEST  SITE,  1971 


Site 

ranges : 

1=0. 

63  km, 

2  =  1.4  km, 

3=2.5  km, 

4  =  5. 

0  km, 

5  =  7.5  km, 

6  =  10.0  km 

RUN 

TIME 

SITE 

RUN 

TIME 

SITE 

NO. 

EATE 

(PDT) 

NO. 

NO. 

DATE  (PDT) 

NO. 

46 

7/14 

0345 

4 

76 

7/15  1638 

6 

47 

0403 

3 

77 

1648 

6 

48 

0422 

2 

78 

1657 

6 

49 

0453 

6 

79 

1720 

4 

50 

0500 

6 

80 

1727 

4 

51 

0510 

6 

81 

1751 

2 

52 

0530 

4 

82 

1800 

2 

53 

0546 

4 

83 

7/16  0944 

6 

54 

0556 

4 

84 

1001 

5 

55 

7/15 

0643 

6 

85 

1009 

5 

56 

0709 

6 

86 

1027 

4 

57 

0726 

5 

87 

1043 

3 

58 

0737 

5 

88 

1101 

2 

59 

0757 

4 

89 

1139 

1 

60 

0812 

4 

90 

1304 

5 

61 

0829 

3 

91 

1322 

4 

62 

0847 

2 

92 

1330 

4 

63 

0917 

1 

93 

1355 

3 

64 

0924 

1 

94 

1415 

2 

65 

1008 

6 

95 

1440 

1 

66 

1028 

5 

96 

1448 

1 

67 

1047 

4 

97 

1510 

2 

68 

4 

98 

7/19 

4 

69 

1115 

3 

99 

4 

70 

1132 

2 

100 

4 

71 

1153 

1 

101 

1940 

4 

72 

1509 

2 

102 

1954 

4 

73 

. 

1530 

3 

103 

2011 

4 

74 

1602 

4 

104 

2031 

4 

75 

1620 

5 

105 

2044 

4 
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Table  2  (Concluded) 


DATA  SUMMARY— RECIPROCAL  SLANT  PATH  EXPERIMENT, 
NEVADA  TEST  SITE,  1971 


Site 

ranges 

:  1  - 

0.63  km, 

2  ■  1,4  km, 

3  -  2 

,5  km, 

4  ■ 

5,0  km, 

5  ■  7.5  km, 

6  ■  10 

.0  km 

RUN 

TIME 

SITE 

RUN 

TIME 

SITE 

NO. 

DATE 

Ct>DT) 

NO, 

NO. 

DATE 

(PDT) 

NO. 

106 

7/19 

2106 

4 

136 

7/21 

1156 

4 

107 

2124 

4 

137 

1204 

4 

108 

2195 

6 

138 

1226 

5 

109 

2204 

6 

138A 

1234 

5 

110 

2227 

5 

139 

1255 

6 

140 

1408 

1 

111 

2247 

4 

112 

2257 

4 

141 

7/22 

0900 

2 

113 

2308 

4 

142 

0910 

2 

114 

2331 

3 

143 

0930 

3 

115 

2303 

2 

144 

0957 

4 

145 

1014 

n 

116 

0019 

1 

117 

0028 

1 

146 

1036 

6 

118 

0133 

6 

147 

1101 

3 

119 

0156 

5 

148 

1119 

4 

120 

0204 

5 

149 

1134 

5 

150 

1141 

5 

121 

0222 

4 

122 

0241 

3 

151 

1159 

6 

123 

0250 

3 

151A 

1208 

6 

124 

0312 

2 

152 

1247 

3 

125 

0319 

2 

153 

1256 

3 

154 

1314 

4 

126 

0344 

1 

154A 

1321 

4 

127 

0445 

6 

155 

1337 

5 

128 

0457 

6 

129 

0518 

5 

156 

1354 

6 

130 

0540 

4 

157 

1401 

6 

158 

1409 

6 

131 

0548 

4 

159 

1434 

3 

132 

7/20 

0007 

3 

160 

1443 

3 

133 

7/21 

1029 

1 

134 

1102 

2 

161 

1500 

4 

135 

1132 

3 

162 

1517 

5 

163 

1534 

8 

164 

1543 

6 
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In  praotioa,  tha  laaer  la  aligned  with  tha  raoaivar  before  tha 
atari  of  tha  run.  Tha  run  than  bafina  with  tha  reoording  of  calibration 
lava  la,  aftar  which  tha  laaar  output  ia  ahielded  for  approximately  IS  a 


while  tha  noiaa  in  tha  raoaivar  ayatem  ia  raoordad.  Tha  PDA  ia  than 


stoppad,  and  ita  output  ia  diapiayad  on  tha  x-y  recorder  while  tha  data 
oontinua  to  ba  raoordad  on  tape.  Next,  tha  PDA  oollaota  tha  danaity  of 
tha  aointillation  data,  and  thia  ia  plotted  over  tha  calibration  plot. 

An  example  of  a  PDK  ooileoted  on  line  ia  Bivan  in  figure  11.  Tha 
figure  illuatrataa  tha  calibration  and  "noiaa"  oheoka  made  at  tha 
beginning  of  tha  run,  in  addition  to  tha  aignal  reoeived  along  a  ground- 
to-air  propagation  path. 


IV  DATA  REDUCTION  AND  ANALYSIS 


A.  Data  Proofing 

Aa  dlaouiaed  in  Chapter  III,  th«  laair  aointi nation  and  temperature 
fluctuation  data  wara  raoordad  in  analog  faahion  on  FW  magnatio  tapa. 

Tha  probability  danaity  analyaar  (PDA)  diaouaaad  aarliar  waa  uaad  to 
prooaaa  tha  data  to  obtain  frequency  distribution* |  tha  varianoa  of  tha 
signal  la  than  aaslly  obtalnad  for  una  in  tha  application  of  tha  thaorlaa 
of  Tataraki  and  Kolmogoroff  (1041), 

In  oparation,  tha  PDA  takas  a  Urge  numbar  of  samples  of  tho  analog 
wavaform  and  raoords  tha  numbar  of  aamplaa  that  ooour  in  aaoh  of  1024 
voltage  bins.  Typloally,  about  ona  million  aamplaa  are  obtalnad,  Tha 
data  are  presented  aa  a  relative  frequency-of-occurrence  histogram  that 
constitutes  an  estimate  or  "measured  value"  of  tho  probability  danaity 
function  of  tha  scintillation  or  thermal  waveform  with  apeotral  content 
from  da  to  12  OH*.  (Actually  tha  spectral  content  of  the  system  Is 
limited  by  the  frequency  passband  of  tha  tapa  reoordor  that  is  about 
0  to  12&0  Ha  at  tha  tap*  speed  used.)  Specifically,  a  sampling  oscil¬ 
loscope  in  th*  PDA  takes  nearly  instantaneous  samples  of  the  signal  and 
generates  a  boxcar  wavaform  whoso  height  la  proportional  to  th*  sampled 
voltage.  In  the  pulse-height  analyser  section,  an  analog-to-dlgitnl 
convertor  categorises  the  section  of  the  waveform  corresponding  to  each 
sample  aa  being  within  one  of  th*  1024  contiguous  voltage  intervals.  The 
number  of  times  a  sample  occurs  in  each  bin  during  a  given  test  is  recorded 
in  memory.  After  a  selected  number  of  samples  has  been  obtained,  the  data 
in  the  memory,  are  displayed  on  a  CRT  and  can  be  analysed  to  obtain  the 
atatistios  of  the  probability  function. 


Precadlng  page  blink 
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For  this  purpose ,  ws  hsvt  constructed  s  hardware  interface  to 

* 

transfer  the  data  to  a  Wang  780/702  Programmable  Calculator/Plotter. 

When  a  density  of  sufficient  length  has  been  collected)  it  is  trans¬ 
ferred  through  the  interface  to  the  calculator.  Software  programs  have 
been  developed  to  analyse  the  density  to  obtain  the  mean,  standard 
deviation*  skewness,  and  kurtosis  of  the  distribution.  Theso  output 
data  oan  then  be  displayed  with  an  on-line  typewriter-plotter. 

Figure  18  is  an  example  of  a  typical  display.  The  probability 
density  histogram  generated  by  the  PDA  is  indioated  by  the  oroaaes  in 
the  figure.  To  decrease  the  time  required  to  print  the  curves,  only 
every  fifth  point  is  plotted j  however,  this  spacing  can  be  controlled  by 
the  operator  according  to  the  resolution  required.  The  distribution  has 
been  normalisod  about  a  aero  mean,  and  the  aeoond,  third,  and  fourth 
momenta  aro  Hated  benoath  the  plot.  The  Gaussian  curve  corresponding 
to  the  variance  of  the  measured  distribution  is  also  plotted  routinely; 
it  is  indicated  by  the  heavy  dots  in  the  figure.  The  difference  between 
the  two  curves  is  shown  by  the  distribution  of  stars  at  the  bottom  of 
the  plot. 

A  powerful  capability  of  the  data  reduction  program  is  illustrated 
in  Figure  13,  At  the  longer  ranges  the  negative,  normalized  intensity 
fluctuations  often  drop  into  the  background  noise  level  of  the  PMT 
(photomultiplier  tube)  detector,  Fortunutoly,  huwover,  most  of  the 
scintillations  are  detected,  although  the  distribution  is  severely  skewed 
in  the  positive  sense.  One  of  the  results  of  our  eurlier  programs 
(.Johnson  ot  ul.,  1970)  was  the  recognition  that  the  scintillation  can  be 
accurately  represented  by  a  log-normal  distribution.  Therefore,  we  have 
developed  un  optional  software  subprogram  thut  automatically  locates  the 

"  Wang  Laboratories 
838  North  Street 
Tewksbury,  Massachusetts 
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PROBABILITY  DISTRIBUTION  WITH  OAUSSIAN  PIT 


FIGURE  12 


PROBABILITY  DISTRIBUTION  WITH  GAUSSIAN 
FIT,  RUN  NO.  "M" 
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FIT,  RUN  NO.  38 
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weighted  median  of  the  distribution  and  computes  the  standard  deviation 
from  the  positive  half  of  the  data.  The  results  of  this  x'outine  are 
illustrated  in  Figure  13  for  an  evening  run  at  a  10-km  range  along  an 
air-to-surface  path.  The  signal  for  this  run  begins  to  drop  off  system¬ 
atically  below  approximately  -3.5  dB.  The  weighted  median  is  computed 
for  the  -3.5-  to  +3.5-dB  range  and  normalized,  or  set  to  zero;  the 
standard  deviation  is  computed  for  the  half  of  the  distribution  where 
the  relative  log  intensity  is  positive.  The  fitted  curve  is  seen  to 
represent  the  "good"  data  accurately. 

During  the  course  of  the  data  processing,  we  also  sought  to  examine 
the  short  term  variations  of  the  variance  of  the  scintillation.  This  is 
especially  important  as  the  scintillation  is  a  result  of  turbulence  con¬ 
ditions  over  the  entire  path,  while  the  temperature  fluctuations  are 
measured  at  discrete  locations.  Figures  14  and  15  illustrate  the  nature 
of  the  scintillation  distribution  measured  along  a  2.5-km  surface-to-air 
path  over  three  consecutive  1-min.  periods  during  inversion  (midnight)  and 
lapse  (late  morning)  conditions.  The  scintillation  for  the  nighttime  runs 
(Figure  14)  is  quite  consistent,  ranging  from  1,7  to  1.9  dB.  The  charac¬ 
teristic  meteorological  feature  was  a  surface-based  radiation  inversion 
o 

of  2.5  C  between  the  surface  and  90  m;  winds  were  light  (1.5  m/s)  although 
steady,  both  in  speed  and  direction,  from  the  SE.  An  interesting  feature 
is  illustrated  by  the  daytime  sequence  (Figure  15):  The  standard  deviation 
of  the  scintillation  intensity  is  quite  steady  about  1.9  dB,  but  a  shift  in 
the  mean  level  of  the  scintillation  intensity  occurs  during  the  second 
minute  [Figure  15(b)],  The  distribution  appears  to  be  a  composite  of  two 
populations  whose  mean  levels  differ  by  about  5  dB.  The  lower  atmosphere 
was  unstable,  with  a  temperature  decrease  of  3.8°C  between  3  and  16  m; 
winds  were  light  (2  m/s)  and  variable  from  the  south.  Moderately  strong 
convective  activity  was  present  under  these  conditions  and  it  can  be 
expected  that  the  horizontal  temperature  structure  of  the  atmosphere  was 
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PROBABILITY  DISTRIBUTION  WITH  GAUSSIAN  FIT 
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FIGURE  14  PROBABILITY  DISTRIBUTION  WITH  GAUSSIAN 
FIT,  RUN  NO.  33 
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Table  3 


RESULTS  FROM  NTS  PATH  GEOMETRY  AND 
RANGE  DEPENDENCE  EXPERIMENT 


T  (®C> 
I  *  3m 


SURFACE 

RECEIVER 

TOWER 

RECEIVER 

1.65 

1.22 

1.34 

0.80 

1.41 

0.97 

2.25 

1.77 

2.88 

3.69 

4.01 

5.30 

4.32 

4.04 

4.72 

5.68 

1.51 

1.15 

2.25 

1.88 

3.02 

2.79 

2.12 

3.24 

2.68 

4.13 

1.52 

1.51 

1.46 

1.74 

1.93 

2.56 

3.22 

4.01 

3.38 

5.26 

2.06 

2.85 

1.69 

2.26 

1.16 

2.00 

0.84 

0.99 

0.25 

0.45 

0.91 

0.87 

1.08 

1.06 

1.05 

1.36 

1.46 

1.62 

.  2.61 

2.66 

2.92 

2.94 

3.79 

3.07 

a • 

(°C) 

x  *  2m 

z  =  92m 

x  =  213m 

x*460m 

0.786 

0.561 

0.596 

0.682 

0.697 

0.608 

0.389 

0.320 

0.465 

0.353 

0.384 

0.246 

0.234 

0.129 

0.082 

0.039 

0.046 

0,117 

0.093 

0.118 

0.162 

0.211 

0.239 

0.2C9 

0.315 


46 
53 
57 
0.034 
0.024 

0.061 

0.047 

0.064 

0.059 

0.072 

0.053 

0.039 

0.028 

0.025 

0.025 

0.025 
.036 
.031 
.022 
.024 


0.040 

0.034 

0.044 

0.024 

0.030 

0.039 

0.057 

0.039 

0.024 

0.031 

0.019 

0.028 

0.018 


0.014 

0.014 

0.012 

0.017 

0.014 

0.015 

0.016 

0.017 

0.014 

0.014 

0.014 
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RESULTS  FROM  NTS  PATH  GEOMETRY  AND 
RANGE  DEPENDENCE  EXPERIMENT 


RUN 

NO. 

\mmmm 

<T 

(°C) 

SURFACE 

RECEIVED 

TOWER 

RECEIVER 

T  («C  ) 
z »  3m 

z :  2m 

z  s  92m 

z  s  213m 

32 

3 

1.47 

1.57 

28.9 

0.152 

0.009 

0.020 

33 

3 

1.98 

2.07 

29.0 

0.116 

0.010 

0.019 

34 

2 

0.63 

0.61 

28.6 

0.174 

KEH3 

0.022 

35 

1 

0.71 

0.48 

28.0 

0.242 

0.016 

0.018 

38 

6 

6.00 

5.88 

26.5 

0.091 

0.069 

0.016 

39 

5 

6.36 

6.53 

26.5 

0.119 

0.029 

0.061 

40 

4 

3.69 

■p|.- 

25.3 

0.099 

0.035 

0.064 

41 

3 

2.38 

0.122 

0.048 

0.054  | 

42 

2 

1.48 

Killit 

25.3 

0.190 

0.028 

0.043 

43 

1 

0.98 

0.63 

26.5 

0.172 

!  0.028 

| 

0.048 

44 

6 

6.24 

5.64 

25.6 

0.165 

|  0.043 

0.058 

45 

5 

5.71 

5.39 

25.0 

0.180 

0.033 

0.046 

46 

4 

4.60 

5.00 

25.3 

0.164 

0.051 

0.042 

47 

3 

2.68 

2.47 

25.5 

0.140 

0.027 

0.042  | 

48 

2 

0.94 

1.05 

24.1 

0.131 

0.025 

0.028  j 

49 

6 

5.56 

8.09 

22.4 

0.288 

0.022 

! 

0.038  1 

50 

6 

6.44 

6.96 

21.7 

0.238 

0.020 

0.034  j 

51 

6 

5.89 

5.62 

22.0 

0.222 

0.054 

0.022  1 

52 

mm 

2.80 

3.03 

21.2 

0.206 

0.043 

0.035 

53 

B 

4.60 

4.94 

22.2 

0.179 

0.061 

0.030 

54 

■ 

4.74 

5.30 

21.9 

0.170 

0.048 

0.037 

55 

6 

4.78 

6.41 

24.1 

0.165 

0.028 

0.027 

57 

n 

1.69 

2.11 

30.3 

0.161 

0.016 

0.017 

58 

mm 

2.40 

2.37 

30.9 

0.166 

0.020 

0.022 

59 

H 

1.40 

2.72 

31.4 

0.273 

0.021 

0.021 

60 

H 

2.01 

2.42 

31.4 

0.325 

0.030 

0.033  ( 

81 

3 

1.37 

1.45 

32.3 

0.236 

0.038 

0.039  ( 

62 

2 

2.06 

1.94 

33.3 

0.386 

0.049 

0.051  C 

63 

Bfl 

•  1.04 

0.92 

34.0 

0.478 

0.064 

0,045  C 

64 

H 

0.88 

34.0 

0.513 

0.052 

0.030  C 

0,045 

0.030 
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RESULTS  FROM  NTS  PATH  GEOMETRY  AND 
RANGE  DEPENDENCE  EXPERIMENT 


RUN 

NO. 

SITE 

NO. 

cr  (4B)* 

a 

(°C) 

SURFACE 

RECEIVER 

TOWER 

RECEIVER 

T  (°C  ) 
z  a  3m 

z :  2m 

z  =  92m 

z  =  213m 

z  -  460  m 

65 

6 

6.07 

6.90 

34.4 

0.579 

0.064 

0.038 

0.069 

67 

4 

2.27 

2.31 

34.2 

0.225 

0.028 

0.025 

0.040 

69 

3 

2.01 

1.78 

36.4 

0.340 

0.083 

0.071 

0.026 

70 

2 

2.37 

2.98 

38.3 

0.588 

0.063 

0.024 

0.037 

71 

1 

0.69 

1.09 

38.3 

0.051 

0.028 

0.043 

72 

2 

1.01 

0.88 

38.3 

0.053 

0.036 

0.033 

73 

3 

2.06 

2.26 

37.1 

0.024 

0.018 

0.019 

74 

4 

2.89 

3.13 

38.3 

BBS 

0.050 

0.040 

M_ 

75 

5 

3.16 

3.65 

35.6 

IBB 

0.011 

0.017 

0.055 

78 

6 

4.08 

4.01 

36.2 

0.268 

0.027 

0.032 

0.015 

79 

4 

2.53 

2.10 

37.9 

0.419 

0.060 

0.031 

0.027 

80 

4 

2.92 

2.17 

37.2 

0.396 

0.049 

0.031 

0.021 

82 

2 

1.63 

1.12 

37.1 

0.435 

0.038 

0.037 

0.026 

85 

5 

4.50 

5.17 

34.2 

0.481 

0.048 

0.044 

0.035 

86 

4 

4.89 

4.88 

33.3 

0.726 

0.075 

0.053 

0.033 

87 

3 

2.77 

3.24 

33.4 

0.570 

0.039 

0.079 

0.022 

88 

2 

1.99 

1.69 

34.4 

0.609 

0.055 

0.084 

0.016 

89 

1 

0.97 

1.12 

34.5 

0,742 

0.035 

0.090 

0.071 

93 

3 

3.57 

4.24 

38.1 

0.489 

0.102 

0.094 

0.082 

94 

2 

1.35 

1.18 

37.5 

0.518 

0.064 

0.075 

0.020 

96 

S  I 

1.40 

1.25 

38.0 

0.809 

0.061 

0.085 

0.045 

97 

1 

2.31 

1.87 

37.8 

0.701 

0.054 

0.063 

— 

101 

■ 

0.58 

0,74 

33.0 

0.045 

0.010 

— 

0.014 

102 

Wmm 

0.42 

0.48 

32.3 

0.044 

0.010 

— 

0.015 

103 

0.48 

0.55 

31.3 

0.057 

0.011 

" 

0.015 

104 

0.39 

0.50 

30.3 

0.087 

0.010 

MM 

0.015 

105 

km 

0.76 

0.64 

29.8 

0.082 

0.011 

MM 

0.015 

106 

pLlI 

1.83 

1.04 

28.7 

0.155 

0.024 

— 

0.013 

107 

WjM 

1.14 

1.18 

28.3 

0.120 

0.019 

MM 

0.014 

108 

6 

2.40 

2.69 

28.3 

0.139 

0.015 

— 

0.014 
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Tabl«  a  (Conoiudad) 


RESULTS  FROM  NT*  RATH  GEOMETRY  AND 
RANOE  DERENDENCE  EXPERIMENT 


j  9  (4  8)  w 

SURFACE 

RECEIVER 

TOWER 

RECEIVER 

4.88 

4.93 

8.88 

6.86 

8.34 

8.67 

4.08 

4.38 

3.84 

4.33 

8.19 

6.37 

5.89 

8.83 

3.60 

3.73 

3.96 

3.87 

6.3$ 

6.88 

6.10 

6.18 

5.48 

6.73 

5.05 

_ 1 

7.89 

T  (•C) 
i  *  Sm 

34.4 

31.8 

33.6 

33.6 

37.3 

36.7 

36.4 

37.1 

37.1 

37.1 


0.786 

0,784 

0.888 

0.608 

0.613 

0.878 

0.668 

0.783 

0.883 

0.788 

0.684 

0.688 

0.733 


9  t»C) 

i»»lm  i  •  Him  i  •  460m 

0.083  0,064  0,073 


0.083 

0.083 

0.107 

0.068 

0.083 

0.071 

0.030 

0.078 

0.088 

0.078 

0.088 

0.114 

0.036 


0.071 


0.048 

0.030 

0.083 

0.048 

0.088 

0.081 

0.108 

0.034 


0,078 

0,047 

0.073 

0.080 

0,080 

0.030 

0.088 

0.088 

0.043 

0.034 

0.038 

0.068 

0.030 


*  Divide  values  by  4.34  to  oonvert  from  dB  to  conventional  natural 
logarithms. 

-  Data  miaaing  or  invalid. 


q«Ue  inhomogeneuua,  These  lateral  thermal  gradients  oeuld  cause  the 
Imr  beam  to  wander,  thereby  producing  tha  blmodal  distribution  ahown 
in  tha  figure,  Pigments  of  data  containing  than  blmodal  distributions 
ha va  bean  deleted  in  tha  prooaaaing  routine  by  tha  operator, 

Similarly,  caaaa  of  poor  reoeivtr-tranamitttr  alignment,  iaaar  Mal¬ 
function..  aaohanioal  vibration,  and  ao  forth  hava  baan  detected  and 
umlttad  from  tha  analyaaa«  furthermore,  wa  hava  ohoaan  to  diaoard  thoaa 
runs  whara  only  &  single  Iaaar  may  hava  baan  oparating  aatiafaotorily  or 
where  tha  nanr»aurfaoe  differential  thermometer  waa  inoperative,  Of 
the  total  of  US  rune  made  in  tha  field,  199  wara  suitable  for  aubaequant 
analyaia  aa  dictated  by  tha  above  criteria,  Tha  atandard  davlationa  of 
the  laaer  and  differential  temperature  data  for  theaa  rune  are  tabulated 
in  Table  9, 

B,  Heaulta 

1.  general 

The  analyaia  procedure  deaoribed  in  tha  preceding  aeotlon  pro- 
vldea  valuea  of  the  atandard  deviation  O’)  of  the  log-lntenaity  algnala 
from  the  two  laaer  recolvera  and  of  the  differential  temperature  fluu- 
tuatlona  (AT)  from  the  four  differential  thermometer# • 

The  menaurementa  of  o(AT)  are  uaed  to  derive  valuea  of  C  ,  t he 

n 

rofraotive  index  atruoture  conatant.  Specifically,  the  temperature 

atruoturo  function,  l)„,  ia  defined  aa 

T 

Dt  b  [Td'j)  -  T(ra)  )*  -  (At)  .  (8) 

In  the  inertial  aubrange  of  the  apeotrum  of  atmoapheric  turbulence,  the 
turbulenoe  ia  iaotropic  and  the  Kolmogorov  (1941)  "two-thirda"  law  can 
be  applied  to  describe  the  refractive  index  atruoture  funotion,  D^t 
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D  <r)  ■  C«(r,)  “  n(ra)3*  ■  C  *  rl/a 
n  i  ■  n 


(•) 


where  r  la  tha  separation  distance  and  C  la  the  retr active  Index  atruo- 

n 

tura  aonatanti  Whan  tha  Kolmogorov  law  la  applied  to  tha  temperature 
fluctuations,  tha  analog  of  Iquatlon  (•)  la  aeon  to  ha 


„  „  *  a/a 

0_  ■  C_  r 
T  T 


(10) 


whara  la  tha  temperature  atruo tura  constant,  On  oomblnlng  Equation 
(10)  with  (I), 

c  «  u(6T)r"l/a  ,  (11) 

T 

It  can  aaaily  ba  ahown  (e«g»,  Tatarskl,  1071)  that  tha  atruatura 
oonatanta  for  tha  rafraotlva  Indax  and  tamparatura  ara  aaaantlally  propor¬ 
tional  and i  whan  neglecting  tha  minimal  affaota  of  wator  vapor,  ara 
ralatad  by 

c  ■  (?•  y  10"«)  (p/T®)  CT  ,  (18) 

whara  p  la  atmoaphario  praaaura  In  millibars,  and  T  la  tha  abaoluta  air 
tamparatura  in  degrees  Kalvin,  Combining  Equation*  (11)  and  (18), 


C  ■  (79  X  10-0)  (p/r^r-l/S  o(6 T)  .  (13) 

n 


Surface  praaaura  obaarvatlona  wara  obtainad  from  the  naarby 
National  Weather  Sarvioa  atatlon  at  Yuooa,  Nevada  (elevation  1800  m  mal) , 
Baoauaa  of  tha  large  variation  in  height  over  which  tha  differential 
tamparatura  moaauremant*  ware  made,  it  waa  neceaaary  to  compute  praaaura 
valuaa  at  the  varioua  levels.  in  thla  regard,  the  preaeure-height 
relationship  in  a  standard  atmoaphara  (Barry  at  al.,  1945)  waa  uaad, 
where 
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» 


(14) 


P  -  PQ 


«/K  y 


and  a  la  tha  elevation  above  the  reference  laval  where  temperature  and 

pressure  ara  T  and  p  ,  respectively}  g  la  tha  acceleration  of  gravity | 
oo 

V  la  an  avaraga  tamparatura  lapaa  rata  (assumed  aqual  to  8,9°K  km"1) j 
and  K#  la  tha  gaa  constant.  Undar  tha  aaaumad  conditions,  tha  exponent 
haa  tha  valua  of  9.96. 


For  a  alant  path  (aubaorlpt  a)  of  langth  It  with  tha  ooordlnata 
origin  at  tha  reoeiver,  Tatarakl  (1981)  haa  davalopad  a  thaoratloal 
axpraaalon  for  tha  atandard  deviation  (c^)  of  tha  log  of  tha  Intanalty 
fluotuatlona  of  a  plana  monoohromatlo  baam  In  a  turbulant  atmoaphara, 
whara 

it 

-  1.4.  (Ml ,na[J  c/wr8'*  dr 

O 

For  a  apharloal  wave,  tha  aquation  for  *s  haa  tha  aama  form  but  tha 
constant  haa  tha  valua  0.97. 


1/9 


(19) 


To  avaluata  tha  intagral,  C  must  ba  apaoifiad  aa  a  funotion  of 

n 

distance  along  tha  Nlant  path.  Thara  ia  inauffioiant  knowledge  about  the 

vertical  atructuro  of  and  indeed  tha  measurements  made  during  this 

program  represent  a  substantial  contribution  in  this  area.  In  integrating 

equation  (19),  we  have  assumed  that  C  ia  independent  of  distance  in  the 

n 

horizontal  plane.  Vertical  variations  of  C  have  been  treated  in  two 

n 

waya:  (1)  assuming  a  simple  exponential  decay,  and  (9)  on  the  basis  of  a 
linear  variation  between  measurement  levels. 


The  model  assumption  for  is  that  used  in  our  earlier  studies 
(Johnson  et  ul.,  1970,  and  Dabberdt  and  Johnson,  1971), 


C  «  C  exp  (-ka)  ,  (16) 

n  no 
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where  k  ii  t  selectable  parameter  with  unitv  of  m“l  and  C  is  the  near- 

no 

ground,  measured  value.  Whan  tha  proposed  C  profile  ia  aubatitutad 

n 

into  Equation  (15),  and  geometrical  oonaidarationa  Involving  tha  alant- 
path  elevation  angle  (5)  are  taken  into  aooount,  we  obtain  tha  following: 


ta 


1.43  <2nA)7/la  C  A 
no 


i/a 


(17) 


Tha  value  of  the  oonatant  ia  changed  slightly  from  Tatarakl'a  plane-wove 
value  in  keeping  with  our  earlier  format  as  determined  by  the  reaults  of 
horiaontal  propagation  experiments  (Johnson  at  al.,  1970).  The  range 
Integral  (A)  in  Equation  (17)  ie  dependent  on  the  direction  of  propagation 
and  ia  given  by 


Case  I  (surfaoe  laser-o leva ted  receiver): 

R 


A  «  A 


/-B(R-r) 
* 


rS/6dr  . 


(1.) 


Case  II  (elevated  laser-surf aoe  receiver) : 

R 

_  -Br  5/6 

A  -  A „  -I  a  r  dr  , 


-/ 


(19) 


where  B  ■  2k  sin  9.  Equations  (18)  and  (19)  are  then  solved  by  numerical 
integration. 


In  the  second  method,  Equation  (15)  Is  solved  analytically 
through  the  use  of  a  linear  interpolation  of  C n  between  measurement 
levols,  where 


C 

n 


V 


(a^  ^  1  +  a  [a  -  }  . 


(20) 


The  braoketed  term  in  Equation  (15)  is  then  integrated  using  Equation 
(20) ,  where 
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Case  I  (surfuce  luser-elevated  receiver) : 


**■*!■£  Cn ,  4-i /  [l  *  *4-1  <rfl  -  r>  ■ln  9]SrV6  dr'  (21> 

i*l  ^ 


Cnee  II  (elevated  loser-surface  receiver): 


R  >=  R  » 

m 


1T1 

£  C!l  f  ,  [l  -  *1*!  *  <*1  *ln  9>r] 

1  =  1  B 


5/6  _ 
r  dr. 


The  subscript  notation  1  denotes  the  layer  over  which  the  integration  is 
performed  where  r^  =  a  esc  @,  and  is  the  constant  in  Equation  (20) 
for  the  layer  with  lower  bound  z^  .  The  use  of  Equations  (18),  (19),  (21), 
and  (22)  is  discussed  in  subsequent  sections  of  this  report. 


2.  Height  Dependence  of  the  Refractive  Index 
Structure  Function 

Equation  (13)  is  used  to  compute  values  of  the  refractive  index 

structure  function,  C  ,  using  the  mean  and  fluctuation  temperature  measure- 

n 

meits ,  the  computed  profile  of  atmospheric  pressure  [Equation  (14)],  and 

a  60-cm  separation  distance.  The  C  data  were  stratified  according  to  the 

n 

atmospheric  stability  in  the  lowest  92  m: 

Inversion  conditions 

T(2  m)  -  T(92m)  <  0.9°C 
Lapse  conditions 

T(2  m)  -  T(92  m)  >  0.9°C  . 

The  individual  C  profiles  were  then  used  to  evaluate  the  parameter  k 
n 

from  Equation  (16),  where 

k  =  -  —  log  (C  /C  ) 
z  e  n  no 


Three  values  of  k  were  obtained  for  each  profile  and  these  have  been 
stratified  by  the  stability.  The  mean  and  standard  deviation  of  the 
ratio  cn/cno  ant*  the  average  k  value  for  each  category  are  summarized 
in  Table  4.  Only  those  runs  were  used  where  simultaneous  AT  measure¬ 
ments  were  obtained  at  all  four  heights;  these  include  30  inversion  and 
68  lapse  cases. 

Tab le  4 

HEIGHT  AND  STABILITY  DEPENDENCE 

OF  C  /C  AND  k 
n  no 


Layer 

Inversion 

Lapse 

<m) 

<c ~7c  ) 

n  no 

k  [m-1] 

(C  /C  ) 
n  no 

k  [m“l] 

0-92 

0.179  ±  0.14 

0.0187 

0.129  ±  0.08 

0.0223 

0-213 

0.194  ±  0.13 

0.0077 

0.109  ±  0.07 

0.0104 

0-460 

0.106  ±  0.04 

0.0049 

0.086  ±  0.07 

0.0053 

The  observed  height  dependence  of  C  is  illustrated  in  Figure  16, 

n 

together  with  the  theoretical  structure  given  by  Equation  (16)  using  two 

typical  values  of  k.  It  is  highly  significant  that  the  theoretical, 

exponential  decay  of  C  with  height  does  not  provide  a  particularly  good 

n 

representation  of  the  details  of  the  average  profile  structure;  individual 
profiles  are  even  less  well  represented.  The  reasons  for  the  apparent 
failure  of  this  simple  model  to  simulate  observed  conditions  can  be 
explained  through  a  consideration  of  both  micro-  and  mesometeorologlcal 
processes.  In  the  case  of  the  latter,  the  effects  of  such  meteorological 
phenomena  as  frontal  activity,  squall  lines,  and  air  mass  structure  must 
be  considered  when  describing  the  temporal  and  spatial  variations  of  C  . 


FIGURE  18  ASSUMED  EXPONENTIAL  FORM  FOR  VERTICAL 
PROFILE  OF  REFRACTIVE  INDEX  STRUCTURE 
CONSTANT  SHOWING  THE  COMPUTED  AVERAGE 
VALUES  AT  92,  213,  AND  480  m  FOR  BOTH 
INVERSION  (-*-1  AND  LAPSE  ("°H  CONDITIONS 


Clearly,  these  cannot  be  simulated  on  the  basis  of  purely  local  conditions. 

In  the  absence  of  these  effects,  micrometeorological  influences  will 

determine  the  C  structure.  Temporal  variations  will  then  occur  with  a 
n 

predominantly  diurnal  periodicity.  Shortly  after  sunrise,  solar  energy 
will  be  absorbed  at  the  earth's  surface  and  subsequently  partitioned,  in 
part,  as  a  flux  of  sensible  heat  in  both  the  soil  and  atmosphere.  In 
comparison,  direct  heating  of  the  air  by  solar  absorption  is  minimal  and 
can  be  ignored.  Under  these  conditions,  the  input  of  heat  to  the  atmos¬ 
phere  is  controlled  by  the  surface ,  and  the  transport  of  heat  to  higher 
levels  occurs  with  a  finite  time  lag.  Moreover,  the  time  lag  is  height 
dependent.  At  night,  the  situation  is  reversed  because  of  radlational 
cooling  at  the  surface  and  the  subsequent  transport  of  heat  is  from  the 

atmosphere  to  the  surface.  Since  the  magnitude  of  C  is  proportional  to 

n 

the  absolute  value  of  the  atmospheric  sensible  heat  flux,  one  must  con¬ 
sider  both  the  nature  of  diurnal  variations  at  the  surface  and  the  time- 
height  dependence  of  the  atmospheric  diffusion  process. 

As  a  first-order  approximation,  the  problem  may  be  considered  as 
an  analogy  to  the  classical  case  of  thermal  diffusion  of  heat  in  a  homo¬ 
geneous  medium  with  a  periodic  source  function.  Therefore,  the  ne*r- 

surface  time  variation  of  C  is  given  by 

n 

C  (t)  *  Y,  {c  +  AC  cos  <u>  +  0\  ,  (23) 

no  l  no,i  i  no  1  J 

where  the  overbar  denotes  the  time-averaged  value,  i  is  the  number  of  the 
harmonic  period  with  frequency  u>,  C  »  constant  phase  lag,  and  t  is 
time.  Similarly,  the  nonsurface  variation  is  expressed  as 

C  (*,t)  ■  Y*(c  +  AC  cos  (at  -  s/Z,  +  •  (24) 

n  f  1  n,l  in  i  i  J 
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Z  is  the  so-called  "damping  depth,"  where  the  amplitude  of  the  wave  is 
1/e  (about  0.37)  times  ‘•'xe  surface  value.  For  illustration  purposes, 
we  may  consider  only  the  diurnal  period  and  make  the  additional  model 
assumptions : 

C  =  C  exp  (-z/Z),  and  (25) 

n  no 

AC  =  AC  exp  (-z/Z).  (26 

n  no 

Combining  Equations  (23-26) ,  the  ratio  C  /C  is 

n  no 


C 

— —  =  exp  (-z/Z) 
no 


+  AC 

no 

+  AC 

no 


cos  (irt-z/Z  + 


cos  (nit  +  Q 


(27 


Although  the  NTS  experimental  program  did  not  encompass  any 
complete  diurnal  cycles,  we  may  take  advantage  of  the  similar  meteoro¬ 
logical  conditions  that  persisted  over  the  first  few  experimental  days 
to  construct  a  "typical"  diurnal  variation  using  data  from  the  following 
periods:  1000-1930,  12  July;  1900,  13  July  -  0600,  14  July;  and  0040-1000, 

15  July.  The  composite  diurnal  cycle  for  is  illustrated  in  Figure  17. 

The  value  of  C  is  about  2.5  x  10"7  m “1/3  while  AC  is  on  the  order  of 
no  no 

2.0  x  1'0“7,  These  values  and  a  100-m  damping  depth  have  been  used  to 

illustrate  the  simple,  descriptive  model  [Equation  (27)]  in  Figure  18, 

Vortical  C  /C  profiles  are  shown  for  six  times  during  the  period;  the 
•n  no 

phase  lag  £  would  be  about  12  hours  so  that  the  profile  at  t  =  0  would 
correspond  io  120C  LB T, 
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3.  Slant-Path  Scintillation. 


The  theoretical  value  of  the  standard  deviation  of  the  log  of 
the  intensity  fluctuations  was  computed  using  Equation  (15)  with  the  two 
model  assumptions  [Equations  (16)  and  (20)]  for  the  height  variation  of 
the  refractive  index  structure  function: 

r  ^  -i  1/2 

ots  =  1.43  (2tt  /\)7/12\j  (?)  r5/6  drj 

o 

C  =  C  exp  ( -kz) ,  and 
n  no 

C  (z  -*  z  )  =  C  (z  )  /  1  +  a  [z  -  z  ]  1 

n  1  2  nil  1  J 


In  making  the  computations  with  the  linear  interpolation  method 

for  C  ,  values  of  the  coefficient  a  were  determined  for  each  run  between 
n  — 

2  and  92  m,  92  and  213  m,  and  213  and  460  m,  and  Equation  (15)  was  then 
integrated  anaytically  using  Equations  (21)  and  (22).  Figures  19  and  20 
depict  the  relationship  between  the  theoretical  and  measured  values  of 
a  along  air-to-ground  and  ground-to-air  propagation  paths,  respectively. 
These  data  have  been  fitted  by  a  nonlinear,  least-squares  computer  program 
to  an  empirical  equation  of  the  form 


a  = 


m  "  1  +  aa  0 

v 


(28) 


where  the  subscript  m  refers  to  the  measured  value,  and  cx  and  3  are  con¬ 
stants.  Equation  (28)  has  been  evaluated  independent  of  the  atmospheric 
stability  and  is  plotted  along  with  the  data  in  the  two  figures.  For  the 
air-to-ground  path,  o.  =  0.47  and  0  •  1.10,  while  for  the  ground-to-air 
case,  &  ■  2.31  and  0  =  0.53. 
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FIGURE  19  SCINTILLATION  ALONG  AN  AIR-TO-GROUND 
PATH  FOR  LAPSE  (+>  AND  INVERSION  (A) 
CASES;  LINEAR  INTERPOLATION  MODEL  FOR 
Cn  (z) —EQUATION  20 


FIGURE  20  SCINTILLATION  ALONG  A  GROUND-TO-AIR 
PATH  FOR  LAPSE  (4)  AND  INVERSION  (A) 
CASES;  LINEAR  INTERPOLATION  MODEL  FOR 
C„  (*) — EQUATION  20 
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Figures  19  and  20  also  illustrate  the  apparent  dependence  of  the 

scintillation  on  the  class  of  atmospheric  stability.  Those  experimental 

runs  that  were  conducted  under  stable  lower  atmospheric  conditions  are 

depicted  by  the  triangular  symbols  in  the  figures;  unstable  conditions 

are  shown  by  the  crosses.  In  interpreting  these  data,  the  theoretical 

sigma  values  provide  a  convenient  scaling  factor  of  atmospheric  (C  )  and 

n 

geometric  <r,3)  conditions.  Clearly,  the  measured  scintillation  is 
consistently  larger  for  stable  conditions  than  for  the  unstable  cases. 

To  examine  these  differences  in  greater  detail,  the  scintillation 
data  have  been  stratified  according  to  both  the  path-type  and  stability 
category  in  Figures  21  through  24.  The  values  of  the  constants  Q,  and  0 
from  the  regression  Equation  (28)  for  each  case  are  summarized  in  Table  5; 
the  individual  functions  are  plotted  on  the  respective  figures. 

Table  5 

VALUES  OF  a  AND  0  FROM  EQUATION  (28) 

USING  THE  MODEL  ASSUMPTION 

C  (z  -  z  )  a  C  (z,)  (l  +  a  [z-z  ]  "l 
n  1  2  nil  1/ 


Atmospheric  Stability 

Path  Type 

Inversion 

Lapse 

All  Cases 

a 

0 

a 

0 

a 

0 

Air-to-ground 

0.16 

1.64 

0.72 

0.82 

0.47 

1.10 

Ground-to-air 

1.38 

0.61 

3.21 

0.42 

2.31 

0.53 

Within  the  broad  range  of  experimental  conditions  represented 
by  these  data,  there  are  indications  of  scintillation  saturation  for 
both  propagation  paths.  While  saturation  may  be  inferred,  the  data 
provide  no  evidence  of  scintillation  superssturation  as  had  been  noted 
earlier  by  Dabberdt  and  Johnson  (1971)  for  near-ground  horizontal  paths. 
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FIGURE  21  SCINTILLATION  ALONG  A  GROUND-TO-AIR  FATH 
FOR  LAF8E  CASES;  LINEAR  INTERPOLATION 
MODEL  FOR  CJi)— EQUATION  20 


tttYHVOMTICAU 


•A-t3«l>30 


If  saturation  dote  indeed  occur  along  th*  slant  path,  than  on* 
would  aapaot  to  find  avidanoa  firat  along  th*  ground-to-air  path  (bacauaa 
of  tha  dominaiict  of  conditions  near  to  th#  •ouroa),  Presumably,  satura¬ 
tion  should  also  firat  Pa  notad  during  unstable  atmoapharln  oonditiona 
whan  tha  magnitude  of  tha  refractive  index  atruotur#  funotion  la  largaat. 
Thaaa  axparimantal  oonditiona  ara  th*  baaia  for  th*  data  presented  in 
Figure  II.  Aa  expected,  thara  appears  to  b*  a  aystamatlo  lava  ling  off 
of  th*  obaarvad  aointiiiation  <*  «*  1.7)  at  larga  dlatanoaa  and  high  C 

W  fi 

valuaa  (aa  raflaotad  in  tha  aoal*  factor ,  u^).  However,  thara  ara  too 
faw  data  pointa  in  thia  ranga  to  fully  oontirw  that  aaturation  haa  lndaad 
oaourrad. 

A  aurpriaing  aapact  of  ground-to-air  aointiiiation  la  iiluatratad 
in  Figure  II.  Thaaa  data  wara  obtalnad  undar  invar* ion  fMghttlma)  con¬ 
dition*  and  consequently  span  a  smaller  rang*  of  o  valuaa  (oa  ahown  in 

% 

Figure  17).  In  aplt*  of  th*  smaller  valuaa  of  a  ,  th*  magnitud*  of  th* 

t 

aQlntlllation  ia  tha  nan*  us  that  obaarvad  undar  th*  unatabl*  atmoaphario 
condition*.  Comparing  tha  two  ossa*,  it  appears  that  th*  aointiiiation 
undar  atabl*  condition*  ia  oomparabl*  to  th*  "unatabl*"  aointiiiation 
whara  acala  valuaa  (j  )  ara  about  two  timaa  largar.  A  aignlficant  con- 

a 

sequence  of  thia  obaarvntion  may  b*  th*  dapandano*  of  th*  aointiiiation 
on  tha  frequency  or  aoal*  of  th*  thermal  fluotuationa. 

Piguraa  13  and  84  ahow  a  similar  pattarn  tor  th*  variation  of 
aointiiiation  along  th*  air-to-ground  path  for  lapa*  and  inversion  con¬ 
ditions,  respectively.  Surprisingly,  thara  is  an  even  stronger  indi¬ 
cation  of  possible  aointiiiation  anturstion  for  the  unstable  casts  with 
this  path  configuration  than  was  notad  earlier  for  tha  ground-to-sir, 
unatabl*  oases.  Again,  th*  small  number  of  oaaev  praoludo  a  definitive 

analysis  of  th*  aaturation  phenomenon.  Th*  maximum  a  value*  are  about 

m 

1.4,  in  contrast  with  th*  ground-to-air  value  of  1.7.  Th*  stability 
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dependence  of  air-to-ground  scintillation  in  i  l  lustra tad  by  the  com¬ 
parison  of  Figure*  83  end  14.  A*  with  the  ground-to-air  ossa*  ■  the 
maximum  observed  eointi  nation  of  about  1.4  for  the  air-to-ground  con¬ 
figuration  ia  Independent  of  stability  (i.e.,  lapae  or  Inversion  con¬ 
ditions).  A*  noted  before,  the  relative  '  dependence  of  a  ia  virtually 

«  m 

identic* 1  for  both  stability  types,  hut  the  comparable  a  values  are 

\ 

again  about  two  times  larger  for  the  lapse  oases.  This  suggests  that 
there  la  a  scaling  effeot  not  only  of  range,  geometry,  wavelength,  and 
turbulenoe  intensity,  but  also  of  the  stability  types,  i.e.,  lapse  or 
inversion.  We  therefore  propose,  in  way  of  conjecture,  that  the  scin¬ 
tillation  magnitude  is  also  dependent  on  the  frequency  (or  scale)  as 
well  as  the  magnitude  of  the  thermal  fluctuations. 

During  lapse  conditions,  there  la  a  decrease  in  the  atmospheric 
potential  temperature  with  height  and  a  tendency  for  conveotive  activity 
resulting  from  the  unstable  density  stratify  cation;  for  inversion  conditions 
the  reverse  is  true  and  the  atmosphere  is  regarded  aa  stable,  Typloally  the 
speatrum  of  near-surf ace  thermal  fluctuations  during  inversion  conditions 
(typioal  of  ,  tear  nights)  ia  characterised  by  relatively  high  frequencies 
(and  low  wavelengths)  due  to  the  height  variation  of  temperature  and  the 
action  of  machanlaal  turbulence.  During  lapse  conditions  (characteristic 
of  sunny  days),  the  spectrum  of  thermal  fluctuations  is  often  augmented 
by  the  effeots  of  buoyency  (free  oonveetlon).  The  impact  of  this  oan- 
vectivs  activity  ia  reflootad  in  the  occurrence  of  a  second  peak  in  the 
spectrum  at  lower  froquenciea  (higher  wavelengths).  In  summary,  nighttime 
condition*  are  characterised  by  relatively  small  though  frequont  turbulent 
elements,  whereas  daytime  conditions  with  free  convection  may  be  dominated 
by  the  larger,  infrequent  convective  cells.  As  a  consequence,  it  may  be 
that  the  "effective  range"  (or  number  of  refractive  occurrences  along  the 
path)  is  significantly  less  for  these  daytime  conditions  although  the  value 
of  may  not  differ  from  the  nighttime  case.  The  idea  of  an  effective 
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rant*  introduced  to  apeolfy  tha  relationship  batwaan  tha  soale  of 

turbulanoa  and  tha  sointl  nation  i  For  equal  valuaa  of  C  and  a  fivan 

n 

propagation  path,  daytime  scintillation  undar  fraa  oonvaotlon  ia 
obaarvad  to  ooour  aa  if  tha  ranga  wara  effectively  raduoad. 

Aa  mantionad  earlier,  tha  measured  aolntlllation  (c  )  rapraaanta 

In 

a  line-integral  valua,  whiia  tha  thaoratioal  valuaa  ara  oonputud  on  tha 

baaia  of  a  faw  point  measurements  of  tha  rafraotiva  indax  atruotura 

funotion,  Two  poaaibla  aouroaa  of  unoartainty  result  t  (1)  tha  vertical 

prof i la  of  C  obtainad  from  tha  maaauramanta  may  not  provtda  auffioiant 
n 

apatial  raaolution  undar  all  atmoaphario  oonditiona,  Thia  la  illuatratad 

by  tha  hypothatloal,  time-dependent  vartloal  profilaa  of  Cr  illuatratad 

in  Figuro  18|  and  (8)  tha  oonoapt  of  "froaen"  turbulanoa  may  not  alwaya 

ba  fulflllad,  l.e,,  tha  time-average  of  C  at  a  point  may  not  alwaya 

n 

provida  an  adaquata  maaaura  of  tha  apaoa  avaraga  of  at  a  givan  tima. 
Whllo  thaaa  affaota  praaumably  add  aoma  noiaa  to  tha  ayatam,  thay  do 
not  affaot  tha  ovarall  or  avaraga  piotura  aa  praaantad  in  Figuraa  81 
through  24, 

At  praaant,  thara  la  lnauffioiant  information  to  paramatariae 

tha  amplituda  and  phaaa  funotlona  in  Equation  (87)  for  tha  Incorporation 

into  tha  eye-safety  avaluation  prooadura  (Chapter  v).  Similarly, 

Equation  (80)  cannot  ba  appliad  in  thia  manner  becauae  it  requirea  in 

aitu  meaaurementa  of  tha  vertical  atruotura  of  C  ,  The  uaa  of  Equation 

A 

(16)  therefore  ia  retained  in  tha  evaluation  procedure  since  it  provldea 
conaarvative  (aafe)  results  and  represents  a  rcodou*ly  successful  simu¬ 
lation  of  conditions  and  is  easily*  parameterized .  Figures  25  through  28 

illustrate  tha  a  -o  relationship  incorporating  Equation  (16)  for  the 
m  t 

data  stratified  in  an  identical  manner  to  tha  presentations  in  Figures 
81  through  84,  rospaotively.  Tha  relative  features  are  the  same  for 
the  two  seta  of  data  and  only  the  absolute  values  of  >j  differ.  The 

X 

unstable,  ground-to-air  cases  (Figure  25)  are  similar  to  the  slant-path 
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FIGURE  26  SCINTILLATION  ALONG  A  GROUND-TO-AIR 
PATH  FOR  INVERSION  CASES;  EXPONENTIAL 
MODEL  FOR  C„(z)— EQUATION  16 
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FIGURE  27  SCINTILLATION  ALONG  AN  AIR-TO-GROUND 
PATH  FOR  LAPSE  CASES;  EXPONENTIAL 
MODEL  FOR  Cn(z) — EQUATION  18 


FIGURE  28  SCINTILLATION  ALONG  AN  AIR-TO-GROUND 
PATH  FOR  INVERSION  CASES;  EXPONENTIAL 
MODEL  FOR  C„<*)— -EQUATION  16 
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data  collected  previously  by  Dabberdt  and  Johnson  (1971);  the  results 
from  the  two  programs  agree  quite  well.  For  the  present  study, 
evaluation  of  the  regression  Equation  (28)  for  these  cases  yields 
values  for  the  two  constants  of  &  =  2.07  and  (3  =  0.45,  while  the 
earlier  (1971)  results  give  &  =  1.79  and  0  =  0,69. 

In  summary,  the  use  of  the  linear  interpolation  model  for 
C^(z)  —  Equation  (20)  ~  is  to  be  preferred  for  the  evaluation  of 
the  field  data  as  it  provides  the  best  representation  of  ambient  con¬ 
ditions  as  determined  from  the  in  situ  measurements.  However,  the 

exponential  model  for  C  (z)  —  Equation  (16)  —  m  >re  suited  for 

n 

incorporation  into  the  eye-hazard  evaluation  procnivra  as  it  is  more 

readily  parameterized.  It  is  desirable,  therefore,  to  compare  the  two 

model  assumptions  in  terms  of  the  o  values  that  result  from  their 

t 

application  in  Tatarski's  theory,  Equation  (8).  Figures  29  through  32 
illustrate  these  a  values  for  both  path  types  and  the  two  stability 
classes.  For  the  two  stability  classes  and  the  ground-to-air  path,  the 
"exponential"  model  consistently  gives  v  values  about  twice  those 

v 

obtained  from  the  "linear"  model,  while  along  the  air-to-ground  path  the 
two  models  provide  essentially  similar  values  of  <7  for  both  stability 
categories,  As  would  be  expected,  the  scatter  is  larger  for  the  inver¬ 
sion  cases  than  with  the  lapse  cases  and  probably  reflects  the  greater 
atmospheric  homogeneity  associated  with  the  latter.  The  application  of 
Equation  (16)  in  the  eye-hazard  evaluation  procedure  has  been  revised 
on  the  basis  of  the  path  dependence  of  the  exponential  model  assumption 

for  C  (z) . 
n 

The  classical  slant-path  propagation  theory  of  Tatarski  [see 

Equation  (1ft)]  predicts  scintillation  values  that  are  heavily  dependent 

on  the  direction  of  propagation  in  the  planetary  boundary  layer  because 

of  the  height  dependence  of  C  .  While  firat  seeking  to  confirm  the 

n 

validity  of  this  relationship,  the  ratio  of  measured  scintillation  along 
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FIGURE  29  GROUND-TO-AIR  PATH— LAPSE  CASES 
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FIGURE  30  GROUND-TO-AIR  PATH— INVERSION  CASES 
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FIGURE  31  AIR-TO-GROUND  PATH— LAPSE  CASES 
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FIGURE  33  AIR-TO-GROUNO  PATH— INVERSION  CASES 
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the  two  reciprocal  paths  was  compared  with  the  theoretical  ratio  using 

Equation  (15)  and  the  C  computations.  A  surprising  result  of  this 

n 

comparison  is  that  there  is  virtually  no  systematic  relationship  that 

can  be  observed.  Therefore,  the  observed  scintillation  ratio  was 

examined  in  view  of  the  corresponding  value  of  the  theoretical,  ground- 

to-air  scintillation  (a  ).  The  use  of  a  is  Justified  on  the  basis  of 

t  t 

its  role  as  a  convenient  scaling  parameter  for  range  and  C^.  The  results 
of  this  comparison  are  illustrated  In  Figure  33  and  have  been  stratified 
according  to  atmospheric  stability.  The  best-fit  linear  regression 
equation  for  each  case  is  also  illustrated.  It  can  be  noted  that  there 
is  a  slight  dependence  of  the  ratio  of  the  magnitude  of  a ,  although  the 
scatter  is  quite  large  and  particularly  so  at  the  low  a  values.  Again, 

v 

the  stability  dependence  is  evident  in  the  difference  in  the  magnitude  of 
the  slope  of  the  regression  line  for  the  two  stability  classes. 

We  may  conclude  that  the  simple  relationship  between  scintilla¬ 
tion  and  path  type,  as  In  the  theory  of  Tatarskl,  is  not  evidenced  by 
our  observations  and  that  additional  research  into  the  complexities  is 
required . 
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FIGURE  33  COMPARISON  OF  THE  RATIO  OF  MEASURED 
VALUES  OF  SCINTILLATION  ALONG  AIR-TO- 
GROUND  AND  GROUND-TO-AIR  PATHS,  AND 
THE  SCALING  FACTOR  (o,|  FOR  DIFFERENT 
STABILITY  CATEGORIES 
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V  LASER  EYE-SAFETY  GUIDELINES 


A .  General 

In  three  earlier  reports  on  this  project  (Johnson  et  al.,  1968, 
1970,  and  Dabberdt  and  Johnson,  1971),  a  procedure  was  described  for 
estimating  the  probability  of  eye  damage  from  exposure  to  a  laser  beam. 
Similar  work  has  been  reported  by  Deitz  (1968,  1969).  The  eye-hazard 
evaluation  procedure  presented  in  this  section  incorporates  minor 
revisions  that  reflect  the  results  of  our  latest  experimental  work. 

The  equations  and  assumptions  underlying  each  part  of  the  eye-hazard 
evaluation  procedure  (as  represented  by  Figures  34  through  41  and  Tables 
6  and  7  will  be  briefly  explained  in  a  subsequent  section. 

It  should  be  emphasized  that  these  eye-safety  guidelines  should  be 
considered  tentative  and  used  with  caution  until  the  new  experimental 
results  that  they  reflect  can  be  confirmed  for  a  variety  of  laser-beam 
and  propagation-path  configurations. 

B.  Eye-Hazard  Evaluation  Procedure 

The  input  parameters  required  for  the  eye-hazard  evaluation  proce¬ 
dure  are  as  follows: 

<p,  latitude  (deg) 

U,  surface  wind  speed  (knots) 

V,  atmospheric  visibility  (miles) 

N,  cloud  cover  (tenths) 

6,  slant-path  elevation  angle  (deg) 
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R,  range  (km) 


P,  laser  peak  power  (W) 

Q,  laser  pulse  rate  (pulses/s) 
laser  wavelength  (p,) 

j3,  laser  beam  divergence  (mrad) 
d,  minimum  beam  diameter  (m) 

T,  exposure  time  (s) 

Laser  type  (CW,  long  pulse,  or  Q-switched) 

Type  of  slant  path  (ground-based  or  airborne  laser) 
Date  and  time. 


Using  the  above  information,  one  can  estimate  the  probability  of 
eye  damage  by  the  method  outlined  below: 

(1)  In  the  case  of  relatively  low-power  lasers,  go  to  Step  (2). 
For  other  lasers,  skip  to  Step  (3). 

(2)  Divide  the  LASER  POWER  P  (W)  by  the  minimum  beam  area 
nd2/4  (m2),  using  the  MINIMUM  BEAM  DIAMETER,  d  (m) ,  to 
compute  the  maximum  power  density  available  for  eye 
damage  (W/m2).  If  this  value  does  not  exceed  the 
appropriate  MAXIMUM  SAFE  POWER  DENSITY  LEVEL  (I  ) 
given  in  Table  7,  then  the  laser  can  be  considered  safe 
and  the  rest  of  the  procedure  disregarded. 

(3)  Use  TIME  OF  YEAR,  HOUR  OF  DAY,  and  LATITUDE,  $ ,  in 
Figure  34  to  find  SOLAR  ELEVATION  ANGLE,  a.  interpolate 
between  graphs  for  dates  between  those  given.  For 
latitudes  in  the  southern  hemisphere,  reverse  the  dates 
on  graphs  (a)  and  (c) .  For  nighttime,  proceed  to  Step 
(5). 

(4)  Use  a  and  CLOUD  COVER,  N  (tenths),  in  Figure  35  to  find 
INSOLATION  CATEGORY,  S  (strong,  moderate,  or  weak). 
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^  m 


(A)  For  daytime,  use  8  am*  SURFACE  WIND,  U  (knots),  in 

Table  «  to  find  NEAR-GROUND  REFRACTIVE  INDEX  STRUCTURE 
constant,  Cn  (NHW"1'3).  For  nighttime,  use  isst  twu 
columns  of  tabli,  Fiven  N, 

00  For  s  near-ground  horizontal  propagation  path,  proceed 
to  Stsp  (S),  For  s  slant  path,  uss  RANGE  k  (km),  and 
SLANT-PATH  I  Lit  VAT  ION  A  NO  LX,  «,  in  Figure  38  to  find 
HANGS  1NTK0HAL,  A  (*-11/8) . 

(7)  For  a  slant  path,  uss  A,  Cn,  LA8IH  WAVSLSNQTH,  X  (u) , 
propagation  direction,  and  ths  atmospheric  stability 

in  Figure  37  to  find  LOO -INTENSITY  STANDARD  DEVIATION,  it 
Assume  LAPSE  CONDITIONS  during  the  daytime  and 
INVERSION  CONDITIONS  at  night.  [Proceed  to  Step  (»)]. 

(8)  For  a  nenr-ground  horisontal  path,  use  R,  C  ,  and  LASER 
WAVELENGTH,  (u) ,  in  Figure  38  to  find  LOG- INTENSITY 
STANDARD  DEVIATION,  J. 

(9)  Use  ATMOSPHERIC  VISIBILITY,  V  (miles),  LASER  BEAM 
DIVERGENCE,  $  (mi lliradinns) ,  and  R  in  Figure  39  to 
find  NORMALIZED  MEAN  POWER  DENSITY,  1/P  (m”a) . 

(10)  Multiply  I/P  by  LA8ER  POWER,  P  (W) ,  to  obtain  MEAN  POWER 
DENSITY,  I  (W  ilT*), 

(11)  Use  X  and  LASER  TYPE  (Q-switohed,  long  pulse,  or  CW)  in 
Table  7  to  obtain  MAXIMUM  SAFE  POWER  DENSITY,  I  „  (W  m“a) . 

HALO  1* 

For  a  CW  laser,  use  the  following  characteristic  pulse 
times:  (1)  3  to  10  ms  for  u  slant  path  involving  un 
aircraft,  and  fur  static  paths  and  U  >  10  knota,  and 
(2)  10  to  500  ms  for  static  paths  and  U  «  10  knota. 

(12)  Divide  I  .  by  I  to  obtain  SAFB-TO-MEAN-POWKR-DENSITY 

safe  - 


(13)  Use  a  and  I  ,  /I  in  Figure  40  to  obtain  INSTANTANEOUS 
safe 

PROBABILITY  OF  EYE  DAMAGE ,  y. 
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( H)  U«*  V  «h<t  KVKNT  KNttJUKMt'Y,  K  In  KUUIH*  <U  ««•, 

«U»nu»tiv»J,y,  in  Kquntlun  tan),  t«»  uUtnin  1 NYKOHATKI) 
KHGftABIUTY  OK  KYK  UAMAUK,  !',  fur  «lvan  KKKMUHK  TIMK* 

T.  Kwr  puUaU  la**!'*,  aat  K  tqual  tu  Uw  l/HKtt  PUUIR  NAYN, 
Q  (|tuUi»/»)t  Kut*  CW  Untti'M,  a«t  K  aqUAl  lu  th«  value* 

It  van  ha  U»w  i 

Mtatie  path 

u  t  io  knot*)  K  *>  ftoo 
u  ■  10  knutai  K  »  1000 

Aii'ui'uft  path}  K  -  aOQO 


I 

I 
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TA-MT1-4C 


FIGURE  35  INSOLATION  CATEGORY  (S)  AS  A  FUNCTION  OF  SOLAR 
ELEVATION  ANGLE  (a)  AND  TENTHS  OF  OPAOUE  CLOUD 
COV2R  INI 

TabU  6 


RELATION  OF  NEAR -GROUND  (2  m)  REFRACTIVE  INDEX  STRUCTURE 

CONSTANT  (C  )  TO  WEATHER  CONDITIONS 

n 


C  (10-6  m“1/3) 
n 

Daytime  Insolation 
Category  (S) 

Nighttime  Conditions 

Surface 

Thin  Overcast 

Wind 

Strong 

Moderate 

Weak 

or  i  5/10 

£  4/10 

(knots) 

Cloudiness 

Cloudiness 

<  3 

0,60 

0,50 

0.40 

0.30 

0.40 

3-6 

0.50 

0.40 

0.25 

0,20 

0.35 

6-10 

0.40 

0.25 

0.15 

0.10 

0.20 

10-12 

0.20 

0.12 

0.05 

0.05 

0.05 

>  12 

0.12 

0.05 

0.05 

0.05 

0.05 
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FIGURE  36  RANGE  INTEGRAL  (A)  AS  A  FUNCTION  OF  RANGE  (R)  AND 
SLANT-PATH  ELEVATION  ANGLE  (6),  COR  THE  TWO  TYPES 
OF  SLANT  PATHS 
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FIGURE  39  NORMALIZED  MEAN  POWER  DENSITY  (i/P)  AS  A  FUNCTION  OF  BEAM  DIVERGENCE  ANGLE  (0>  AND  RANGE  (R). 
FOR  VARIOUS  VALUES  OF  ATMOSPHERIC  VISIBILITY  (V) 
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C.  Paata  for  Kvo-Ha*ard  Evaluation  Procedure 

1,  Figure  34 

The  curve*  in  this  figure  are  baaed  on  the  following 
(Smithsonian  Meteorological  Table*,  1983)! 

Of  ■  ain”l|ain*ainC  +  coa*eo»Coo*t0.863(h-'ia))  |  ,  (89) 

whore  cr  ■  aolar  elevation  angle 
•  ■  latitude 
C  ■  aolar  declination 
h  n  hour  of  day. 

8.  Flguro  38 

The  rotative  inaqlution  atrength  (S).  on  u  acale  from  xoro  to 
ono,  ia  given  by 

S  -  ain  or  '■3U) 

for  a  cloudless,  nonattenuating  atmosphere.  When  the  cloud  cover  (N)  in 
tenths  la  conuidered  und  the  albedo  of  tho  clouds  la  taken  to  be  0.3, 
Equation  (30)  beuomea 

S  ■  Bin  [0(1  -  0.5N/10)]  ,  (3D 

which  is  tho  basis  for  Figure  35.  The  categories  of  S  (strong,  moderate, 
and  weak  insolation)  were  determined  simply  by  dividing  the  range  of  S 
into  three  equal  intervals:  0  to  0.33,  0.33  to  0.87,  and  0.67  to  1.0. 

3,  Table  6 

This  table  is  baaed  principally  upon  empirical  data  from  our 
experiments  and  those  of  other  investigators,  such  as  Wright  and  Schutz 
(1967),  Fried  et  ul.  (1967),  Davis  (1966),  and  Goldstein  et  al.  (1965). 
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Unfortunately,  most  of  the  values  of  reported  by  these  other  workers 
have  not  been  accompanied  by  a  full  aooount  of  the  concurrent  meteorological 
conditions.  Hence  the  dependence  of  on  weather  conditions  as  given  in 
tho  table  is  based  principally  on  our  measurements.  The  values  in  the 
table  must  be  considered  only  as  approximate,  since  some  extrapolation  was 
required,  and  also  since  at  least  one  important  variable,  ground  condition, 
has  been  omitted  for  the  sake  of  simplicity, 

Tho  basic  structure  of  the  table  is  patterned  after  that  given 
by  Paaquill  and  Meade  (see  Slade,  1908)  for  the  estimation  of  atmospheric 
stability  categories.  Since  in  related  to  the  absolute  value  of  the 
stability,  appropriate  valuos  of  were  substituted  for  those  of 
stubility  category, 

4.  Figure  36 

Equations  (18)  and  (19)  were  used  for  the  computations  for  the 
curves  presented  in  this  figure.  The  value  of  tho  constant  k  was  taken 
to  he  0.01  so  the  eye-hazard  entimntos  would  be  on  the  conservative  (safe) 
side,  as  determined  from  tho  BHEN  Tower  experimental  results. 

8.  Figure  37 

This  figure  was  bared  on  Equation  (17)  for  the  slant-path  con¬ 
figuration,  but  we  have  limited  a  to  a  maximum  value  according  to  the 
direction  of  propagation  and  the  atmospheric  stability.  The  curves  are 
thus  derived  from  the  Tatarski  relationship  but  take  into  account  the 
saturation  effeot  indicated  by  experimental  results  (Chapter  IV),  Since 
the  occurrence  of  scintillation  supersaturation  on  the  slant  path  has  not 
yet  been  adequately  confirmed,  we  have  not  made  any  allowance  for  this. 

As  a  result  of  this  omission,  any  resulting  errors  would  give  eye-hazard 
predictions  on  the  safe  side. 
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6 .  Figure  38 

Thia  figure  use*  the  reiulto  of  Dabberdt  and  Johnson  (1971)  — 
Equations  (10)  and  (10) ,  Seotion  I1I-F— and  incorporates  the  super¬ 
saturation  effuot  as  observed  during  that  horiaontal  propagation  exper¬ 
iment.  The  coefficients  of  their  Equation  (10)  usod  in  the  calculation 

are  those  roported  by  Johnson  et  al.  (1970).  These  values  were  ohosen 

* 

btoaui*  they  give  more  ooneervative  (higher)  estimetes  of  a  then  the 

m 

revised  equation.  The  leveling  off  of  the  scintillation  magnitude  in 
the  far-supersaturation  region  is  also  incorporated  in  the  figure;  again, 
a  safe  value  (0.8)  was  chosen. 

7 •  Figure  39 

This  figure  is  derived  principally  from  considerations  of  the 
laser-beam  geometry,  represented  by  the  equation 

4G 

I/P  -  ,  02) 

TTp  R 

-  2 
where  I  *»  moan  power  density  (W/m  ) 

(assumed  to  be  uniform  across  the  beam) 

P  »  laser  output  peak  power  (W) 

P  =  beam  divergence  (milliradians) 

R  *  range  (km) 

C  a  atmospheric  transmission. 

The  transmission  is  given  by 

G  o  o"aR  ,  (33) 
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whoro  a  is  tho  attonuation  coofficiont  (km"1).  To  u  first  approximation, 
tlio  relation  between  tho  atmospheric  visibility  (V,  miles)  and  tho  atten¬ 
uation  coefficient  is  given  by 

V  =  2.43/a  (34) 

for  visiblo  and  ncar-inf x’nred  wavelengths. 

8.  Table  7 

2 

The  maximum  safe  powor  density  levels  (I  ,  W/m  )  given  in 

ShIO 

this  table  are  those  recommended  by  the  U.S,  Air  Force  in  AFM  161-8, 
Change  1  (1971). 


9 .  Figure  40 


The  curves  in  this  figure  were  computed  by  rnnnns  of  the 
following  equations: 


V 


exp(-x  /2)dx 


(35) 


where  V  is  tho  instantaneous  probability  of  eye  damage,  x  is  a  dummy 
variable,  and 

M  =  0/2  +  a"1 In  (i  ^  /I)  ,  (36) 

safe  * 


where  o  is  the  log  intensity  standard  deviation.  The  development  of 
these  equations,  as  detailed  in  the  last  report  (Johnson  et  al.,  1968), 
involves  the  assumption  of  a  log-normal  probability  distribution  of  in¬ 
tensity  fluctuations. 
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10.  Figure  41 


When  one  looks  at  a  laser  boam  continuously  for  T  seconds,  the 
integrated  probability  of  eye  damage  (l*)  over  that  period  will  be  greater 
than  the  instantaneous  probability  of  eye  damage  (Y).  One  simple  way  to 
assess  this  increased  probability  is  to  assume  that  the  intensity  peaks 
occurring  at  n  single  observational  point  represent  separate  events. 

If  (1-Y)  corresponds  to  the  probability  of  no  damage  for  a  single  event 
(i.e.j  the  probability  that  any  one  intensity  peak  will  not  exceed 
Isnfe^  inte*ratecl  probability  of  no  damage  (l-I")  for  a  given 

number  of  occurrences  (N)  in  the  time  period  (T)  is  given  simply  by 
multiplying  the  individual  probabilities: 

i-r  =  (1-Y)N  , 

or  f  =  1  -  (1-Y) 

Since  N  =  FT,  where  F  is  the  frequency  of  events, 

r  »  i  -  «-y>ft  .  <3,) 

This  equation  forms  the  basis  for  Figure  41.  For  y  <<  1.  Equation  (37) 
may  be  rewritten  in  approximate  form  bb 

T  =  yFT  .  (38) 

The  event  frequency  F  may  be  considered  to  be  equivalent  to 

the  pulse  rate  for  a  pulsed  laser,  or  to  the  scintillation  frequency  in 

the  case  of  a  CW  laser.  If  we  assume  the  turbulent  elements  separated 

by  a  mean  correlation  distance  (P  )  to  move  with  the  wind,  then  to  a 

o 

first  approximation  we  have 


F  =  U  /p 
n  o 


> 
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whoro  U  is  the  wind  component  perpendicular  to  the  propagation  path, 
n 

Since  the  variation  in  p  with  range  that  is  predicted  by  theory  dis- 

o 

agrees  with  that  measured;  it  seems  best  for  our  purposes  to  assume  that 
F  depends  only  on  wind  speed  and  beam  sweep  rate.  On  the  basis  of  our 
measurements;  the  following  values  have  been  selected  as  being  representa¬ 
tive  of  the  cases  given: 

Static  paths 

U  s  10  knots:  F  =  500 

U  >  10  knots:  F  =  1000 

Aircraft  slant  paths:  F  =  5000 

These  coefficients  allow  for  the  maximum  scintillation  frequencies 
observed  during  our  earlier  experiments. 

11.  Setting  a  level  of  Acceptable  Risk 

After  the  probability  of  eye  damage  (I")  is  found;  it  becomes 

necessary  to  compare  this  value  with  a  preestablished  value  (r  ) 

safe 

representing  the  highest  probability  Judged  to  be  an  acceptable  risk. 

The  magnitude  of  T  may  change  as  the  nature  of  the  situation  changes, 

Sfl  1 G 

-3  -5 

but  values  in  the  range  of  10  to  10  (one  chance  of  eye  damage  in  1000 
to  100,000  occurrences)  seem  reasonable.  By  way  of  comparison,  the 
probability  of  an  average  American  being  killed  in  an  automobile  accident 
in  the  year  1968  was  one  in  5600,  or  2.8  X  10-4,  while  the  probability  of 
being  injured  in  an  automobile  accident  was  one  in  100,  or  10~2  (World 
Almanac,  1970). 

D.  Examples  of  Use  of  Eye-Hazard  Evaluation  Procedure 

To  illustrate  the  use  of  the  procedure  given  above,  four  examples 
are  presented  in  Table  8.  The  input  parameters  for  these  examples  were 
selected  to  be  as  representative  and  realistic  as  possible  and  to  cover 
a  variety  of  conditions.  Example  2  represents  the  conditions  of  an 
actual  recent  field  experiment. 
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Table  8 


EXAMPLES  OF  LASER  EYE-HAZARD  EVALUATION 


Parameter 

Example 

1  Example  2  Example 

3  Example  4 

A.  SPECIFIED  PARAMETERS 

Date 

21  Dec 

1  May  21  June  21  June 

Time 

2200 

1200  1500 

1500 

(LST) 

<P 

30°  N 

40° N  20° S 

20°  S 

N 

2 

0  6 

6 

(tenths) 

U 

5 

00 

to 

12 

(knots) 

V 

20 

8  20 

20 

(miles) 

Path  Type 

Slant 

Slant  Horiz. 

Horiz. 

Laser  Site 

Airborne 

Ground  - 

- - 

e 

60° 

20°  0° 

0° 

R 

0.4 

5  1 

6 

(km) 

A. 

0.4880 

0.6943  1.06 

1.06 

<U> 

Q 

— 

0.5  10 

10 

(pulses/s) 

8 

0.5 

0.5  0.5 

0.5 

(mrad) 

d 

2 

-  - 

— 

(lO"Jm) 

P 

1 

2.5  x  107  107 

107 

(W) 

T 

1 

0.1  0.1 

0.1 

(s) 

Laser  Type 

cw 

Q-sw.  Q-sw. 

Q-sw. 

B .  DERIVED  PARAMETERS 

Step  No. 

3 

0/ 

0°  60° 

3C>P 

30° 

4 

S 

---  Strong 

Mod. 

Mod. 

5 

Cn 

0.35  0.40 

0.12 

0.12  (10"6nf 

1/3 ) 

6 

A  ] 

..5  K  103  1  x  105 

— 

(mil/*) 

7,8 

g 

0.3  0.65 

0.65 

1.20 

9 

I/P 

30  6  x  10-2 

6.7  9 

x  10-2  (m“2 

) 

10 

l 

30  1.5  x  106 

6.7  x  107  9 

x  105  (W  m“ 

2) 

11 

*safe  °* 

52  x  10Z  3,53  x  106  i 

.77  x  107  1.77  x  107  (w  m" 

2) 

*2  Isafe//I 

1.73  2.36 

0.26 

20 

13 

Y 

10-2  5  x  10-1 

_ 

CO 

1 

o 

rH 

14  J 

F 

5  x  103  0.50 

— 

10  (sec- 

b 

1 

r  > 

10-1  5  x  10-1 

— 

10“3 

15  Safe?* 

No  No 

No 

Yes 

*  -3 

Based  upon  a  selected  value  of  r  =10 

safe 
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Appendix 

EVALUATION  OF  THE  MAXIMUM  RISK 
OF  LASER  EYE  DAMAGE 
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Appendix 


EVALUATION  OF  THE  MAXIMUM  RISK 
OF  LASER  EYE  DAMAGE 

A.  General 

The  results  in  the  main  body  of  this  report  outline  a  detailed 
procedure  for  the  estimation  of  the  probability  of  eye  damage  from 
exposure  to  a  laser  beam.  These  guidelines  are  the  result  of  develop¬ 
ments  and  refinements  that  have  been  made  during  past  several  years 
(Johnson  et  al.,  1968  and  1970;  Deitz,  1968  and  1969;  and  Dabberdt  and 
Johnson,  1971) .  The  procedure  used  in  these  guidelines  incorporates  the 
simulation  of  local  atmospheric  conditions  and  other  factors  in  some 
detail.  While  the  procedure  is  quite  rigorous  and  the  predictions  agree 
well  with  observations,  the  guide  is  not  particularly  well  suited  for 
ready  application  by  the  novice  in  remote  field  situations. 

Recognizing  this  particular  need,  a  simple  guide  has  been  developed 
to  estimate  the  typical,  "worse-case"  level  of  eye  hazard  using  as  input 
only  the  most  readily  available  information: 

•  Total  laser  output  or  energy/pulse 

•  Air  Force  Permissible  Exposure  Levels  (PEL) 

•  Laser  beam  divergence 

•  Eye  pupil  diameter 

e  Laser  wavelength 

e  Laser  pulse  repetition  rate 

e  Clear  atmospheric  conditions 

e  Total  personnel  exposure  time 

e  Given  total  probability  of  super-PEL  exposure 
that  is  acceptable 

e  Approximate  laser  beam  elevation  angle. 


Pracsding  page  blank 
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To  ensure  accurate  estimates  in  a  minimal  amount  of  time,  the  procedure 
incorporates  only  a  few  simple  graphs,  tables,  and  arithmetic  steps. 

B.  Procedure 


The  probability  of  eye  damage  from  exposure  to  a  laser  beam  under 
worst-case  conditions  is  estimated  on  the  basis  of  the  following 
information  on  environmental  and  operating  conditions: 

Range,  R  (km) 

Laser  peak  power,  P  (W) 

Laser  pulse  rate,  Q  (s  *) 

Laser  wavelength,  X  (^) 

Laser  beam  divergence,  B  (mrad) 

Wind  speed,  U  (kts) 

Exposure  time,  T  (s) 

Ground-level  elevation,  z  (km) 

Laser  mode  (CW,  long  pulse,  Q-switched) 

Time  of  day  (day  or  night) 

Typo  of  path  (horizontal  or  slant) 

The  method  is  illustrated  schematically  in  Figure  A-l  and  detail  t 
the  following  step-by-step  procedure: 

Step  1:  Use  GROUND-LEVEL  ELEVATION,  z  (km);  LASER  BEAM  DIVERGENCE, 

B  (mrad);  LASER  WAVELENGTH,  X  (p) ;  and  RANGE,  R  (km);  in 
Figure  A-2  to  determine  the  NORMALIZED  MEAN  POWER  DENSITY, 
i/P  (m-2). 

Step  2:  Multiply  i/P  by_LA5ER  PEAK  POWER,  P  (W) ,  to  obtain  the  MEAN 
POWER  DENSITY,  i  (W  m"2) . 

Step  3:  Use  LASER  WAVELENGTH,  X  (|i)  and  LASER  MODE  in  Table  A-l  to 
obtain  the  MAXIMUM  SAFE  POWER  DENSITY,  IBafe  (W  m"2) . 

For  a  CW  laser,  use  the  following  characteristic  "pulse" 
times: 
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•  2-10  ms  for  a  slant  path  Involving  an  aircraft,  and 
static  paths  with  U  >  10  knots 

•  10-500  ms  for  static  paths  and  U  <  10  knots. 

Divide  I_of_  by  I  to  obtain  the  SAFE-TO-MEAN-POWER  DENSITY 
RATIO,  I*  (n.d.) 

Use  I*  and  the  PATH  TYPE  (horizontal  or  slant)  3  to  obtain 
the  INSTANTANEOUS  PROBABILITY  of  EYE  DAMAGE,  y. 

Ascertain  the  value  of  the  EVENT  FREQUENCY, :F  (s-1).  "For 
pulsed  lasers,  set  F  equal  to  the  LASER  PULSE  RATE,  Q. 

For  CW  lasers,  set  F  equal  to  the  values  given  below, 

•  Static  path 

U  <  10  knots:  F  =  500 

U  >  10,  knots:  F  ■  1000 

e  Aircraft  path:  F  =  5000 

Step  7;  Determine  the  INTEGRATED  PROBABILITY  of  EYE  DAMAGE,  T,  from 
the  INSTANTANEOUS  PROBABILITY  of  EYE  DAMAGE,  yj  the  EVENT 
FREQUENCY,  F;  and  the  EXPOSURE  TIME,  T,  using  Figure  A-4  for 
Short  Exposure  times  (T  e'  0.1  s)  and  the  following  equation 
for  relatively  long  exposure  times  (T  >  1  s) : 

f  ■  y  F  T 


Step  4: 

Step  5: 

Step  6: 


A-5 


FIGURE  A-l  FLOW  CHART  FOR  -WORST-CASE"  USERS  GUIDE  TO  EVALUATE  LASER  EYE  SAFETY  HAZARD 


BEAM  DIVERGENCE  —  mrad  BEAM  DIVERGENCE  —  mrad  BEAM  DIVERGENCE 


T  »  0,1  MHOS 

TAMIMIH 


FIOURI  A -4  INTIQRATID  FROIAIIUTY  OF  IYI  DAMAQI 
(D  At  A  FUNCTION  OF  IVINT  FMOUINCY 
(F|  AND  INSTANTANIOUB  PROBABILITY  OF 
IYI  DAMAQI  hr) 


C,  Summary 


For  convenience,  these  etepe  ere  eummarieed  in  Table  A-li.  Ficure 
A-8  is  derived  principally  lrom  conalderationa  of  the  laaer-beam 
geometry,  reproaentod  by  the  equation, 


i/p  • 


» 


whore 

*»  2 

I  m  mean  power  density  (W/m  ) 

(assumed  to  be  uniform  across  the  beam) 
P  ■  laser  output  peak  power  (W) 

6  ■  beam  divergence  (mllllradlans) 

R  m  range  (km) 

0  m  atmospheric  transmission. 
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Table  A-l 


SUMMARY  OF  EYE -MAX AND  EVALUATION  PROCEDURE 


Step  No, 

Task 

Reference 

1 

Determine  I/P 

Figure  A-8 

V 

Compute  1 

X  *  P  (I/P) 

3 

Determine  XiafC 

Table  A-l 

4 

Compute  X* 

-  W* 

S 

Determine  y 

Figure  A-3 

6 

Aaoertaln  F 

Step  6 

7 

Evaluate  P 

j  Figure  A-A,  or 

1  T  ■  y  F  T 

Tht  atmoapherlc  tranamiaeion  la  taken  from  tha  work  of  Eltarman  (1963) 
for  a  clear  atandard  atmosphere  aa  a  function  of  wavelength  and  height. 
The  maximum  aafe  power  denalty  levela  given  In  Table  A-a  are  thoae 
referenced  by  the  U.8.  Air  Force  (1971).  For  tha  daytime  oaaea, 
the  eye-pupil  diameter  waa  taken  aa  3  mm,  while  a  7-mm  dlam.  waa 
uaod  for  nighttime  caaea. 

The  curve  In  Figure  A-3  waa  computed  by  meana  of  the  following 
equationat  ^ 

V  “  y*exp(-x2/a)dx, 

M 

where  y  la  the  inatantaneoua  probability  of  eye  damage,  x  la  a  dummy 
variable,  and 

M  -  a/a  +  o"lln  » 
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where  a  U  the  log-intensity  standard  deviation.  For  thia  guide,  g  la 
taken  equal  to  the  aaturatlon  (maximum)  value  of  the  log-intensity 
atandard  deviation,  whloh  la  1,1ft  for  a  near-ground  horlaontal  propaga¬ 
tion  path  and  1,85  for  a  alant  path.  Theae  valuea  have  been  determined 
from  an  extensive  aeriea  of  field  experlmenta.  The  development  of 
theae  equatlona,  aa  detailed  in  a  prevloua  report  (Johnaon  et  al.,  1970), 
entaila  the  aaaumption  of  a  log-normal  probability  dlatributlon  of 
lntenalty  fluotuatlona, 

When  one  looka  at  a  laaer  beam  continuously  for  T  aeoonda,  the 
Integrated  probability  of  eye  damage  (O  over  that  period  will  be 
greuter  than  the  inatantaneoua  probability  of  eye  damage  (y).  One 
aimple  way  to  aaaeaa  thia  lnoreaaed  probability  la  to  aaaume  that  the 
lntenalty  peaka  occurring  at  a  aingle  obaervatlonal  point  repreaent 
aeparate  eventa.  If  (1-y)  oorreaponda  to  the  probability  of  no  dumuge 
for  a  aingle  event  (i.e,,  the  probability  that  any  one  lntenalty  peak  will 
not  exoeed  Igufe),  then  the  Integrated  probability  of  no  damage  (1-D 
for  a  given  number  of  ooourrenoea  (N)  in  the  time  period  (T)  ia  given 
aimply  by  multiplying  the  individual  probabilities: 

l-r  -  (i-y)N 

or 

T  ■  1  -  (l-y)N 

or  (since  N  ■  FT,  where  F  is  the  frequency  of  events), 

FT 

T  -  1  -  (1-Y) 

This  equation  forms  the  basis  for  Figure  A-4,  For  y  «  1,  this  equation 
may  be  rewritten  in  approximate  form  as 

T  ■  YFT  , 
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Ths  event  frequency,  i,  may  be  considered  to  he  equivalent  to  the 

pulse  veto  for  •  pulsed  Inter,  or  to  the  •olntillatiou  frequency  for 

n  GW  lnaer.  If  w*  assume  the  turbulent  element a  separated  by  n  mean 

correlation  diatanoe  (p  )  to  move  with  the  wind,  then  (to  a  first 

o 

approximation)  wt  havo 

r  ■  V*. 

where  U  la  the  wind  oomponont  perpendicular  to  the  propagation  path, 
n 

Since  the  variation  in  p  with  range  that  ia  predicted  by  theory 

o 

disagrees  with  that  measured,  it  seems  beat  for  our  purposes  to  assume 
that  P  depends  only  on  wind  speed  and  beam  sweep  rate,  On  the  basis 
of  our  measurements,  the  following  values  have  been  selected  as  being 
representative  of  the  cases  given: 

Static  paths 

U  <  10  knots;  F  ■  800 

U  >  10  knots;  F  -  1000 

Aircraft  slant  paths:  F  ■  5000. 

These  coefficients  allow  for  the  maximum  scintillation  frequencies 
observed  during  our  experiments. 

After  the  probability  of  eye  damage  (F)  is  found,  it  becomes 

necessary  to  compare  this  value  with  a  preestablished  value  (r  .  ) 

safe 

representing  the  highest  probability  Judged  to  be  an  acceptable  risk. 

The  magnitude  of  r  .  may  change  as  the  nature  of  the  situation  changes, 
“*fe  -3  -8 

but  values  in  the  range  of  10  to  10  (one  chanoe  of  eye  damage  in 

1000  to  100,000  occurrences)  seem  reasonable.  By  way  of  comparison, 

the  probability  of  an  average  American  being  killed  in  an  automobile 

-4 

accident  in  the  year  1968  was  one  in  3600,  or  2.8  x  10  ,  while  the 

probability  of  being  injured  in  an  automobile  accident  was  one  in  100, 

-2 

or  10  (World  Almanac,  1970). 

To  illustrate  the  use  of  the  worst-case  procedure  given  in  Section 
B,  four  examples  are  presented  in  Table  A-3. 
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Table  A-3 


EXAMPLES  Of'  WORST-CASE  LASER  EYE-HAZARD  EVALUATION 


Parameter  Example  1 

Example  2 

Example  3 

Example  4 

A.  SPECIFIED  PARAMETERS 

Date 

21  Dec 

1  May 

21  June 

21  June 

Time 

2200 

1200 

1500 

1500 

(LST) 

Lat 

30°  N 

40°  N 

2CP  S 

2tf>S 

U 

5 

8 

12 

12 

(knots) 

Path  Typo  Slant 

Slant 

Horlz. 

Horlz. 

Laser  Site  Airborne 

Ground 

— 

— 

Elev.  Angle  flQP 

20° 

0° 

0° 

R 

0.4 

5 

1 

6 

(km) 

X 

0.4880 

0.6943 

1.06 

1.06 

(p) 

Q 

— 

0.5 

10 

10 

(pulses/s) 

e 

0.5 

0.5 

0.5 

0.5 

(mrad) 

p 

1 

2.5  x  107 

107 

107 

(W) 

T 

1 

0.1 

0.1 

0.1 

(s) 

Laser  Type  CW 

Q-sw. 

Q-sw. 

Q-sw. 

B .  DERIVED  PARAMETERS 

Step  No. 

1 

I/P  5  x  101 

io-i 

5  x  10° 

io-i 

(m~2) 

2 

I  5  X  101 

2,5  x  10fi  5  x  107 

106 

(W  m~2) 

3 

Isafe  0-52  X  102  3.53  x  106  1,77  x  107  1.77  x 

IO7  (W  m"2) 

4 

I*  1.0 

1.4 

0.35 

18 

5 

v  >10-1 

>10"! 

MO-1 

CO 

1 

o 

H 

6 

F  5  x  103 

0.5 

10 

10 

(s  -1) 

7 

r  >5  x  io2  >io-1 

>10  “1 

10-3 

Safe?*  No 

No 

No 

Yes 

* 

Based  on  a  selected  value 


# 
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